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ABSTRACT 

Observations continue to support the interpretation of the anomalous microwave foreground as 
electric dipole radiation from spinning dust grains as proposed by Draine & Lazarian (1998ab). In 
this paper we present a refinement of the original model by improving the treatment of a number of 
physical effects. First, we consider a disk-like grain rotating with angular velocity at an arbitrary angle 
with respect to the grain symmetry axis (i.e., grain wobbling) and derive the rotational damping and 
excitation coefficients arising from infrared emission, plasma-grain interactions and electric dipole 
emission. The angular velocity distribution function and the electric dipole emission spectrum for 
disk-like grains are calculated using the Langevin equation, for cases both with and without fast 
internal relaxation. Our results show that for fast internal relaxation, the peak emissivity of spinning 
dust, compared to earlier studies, increases by a factor of ~ 2 for the Warm Neutral Medium (WNM), 
the Warm Ionized Medium (WIM), the Cold Neutral Medium (CNM) and the Photodissociation 
Region (PDR), and by a factor ~ 4 for Reflection Nebulae (RN). The frequency at the emission 
peak also increases by factors ~1.4 to ^2 for these media. Without internal relaxation, the increase 
of emissivity is comparable, but the emission spectrum is more extended to higher frequency. The 
increased emission results from the non-sphericity of grain shape and from the anisotropy in damping 
and excitation along directions parallel and perpendicular to the grain symmetry axis. Second, we 
provide a detailed numerical study including transient spin-up of grains by single-ion collisions. The 
range of grain size in which single-ion collisions are important is identified. The impulses broaden 
the emission spectrum and increase the peak emissivity for the CNM, WNM and WIM, although the 
increases are not as large as those due to the grain wobbling. In addition, wc present an improved 
treatment of rotational excitation and damping by infrared emission. 

Subject headings: ISM: dust, extinction — ISM: general — galaxies: ISM — infrared: galaxies 



1. INTRODUCTION 

Diffuse Galactic microwave emission in the 10 - 100 
GHz frequency range carries important information on 
the fundamental properties of the interstellar medium, 
but it also interferes with Cosmic Microwave Background 
(CMB) experiments (see Bouchet et al. 1999, Tcgmark 
et al. 2000, Efstathiou 2003). 

It used to be thought that there were only three ma- 
jor components of the diffuse microwave Galactic fore- 
ground: synchrotron emission, free-free radiation from 
plasma (thermal bremsstrahlung) and thermal emission 
from dust. However, in the range of frequency from 10 
to 100 GHz an anomalous microwave foreground which 
was difficult to reconcile with the components above was 
first reported by Kogut et al. (1996a, 1996b). 

de Oliveira-Costa et al. (2002) gave this emission the 
nickname "Foreground X" , reflecting its mysterious na- 
ture. This component is spatially correlated with 100 
/im thermal emission from dust, but its intensity is much 
higher than one would expect by extrapolating the ther- 
mal dust emission spectrum to the microwave range. 
Draine & Lazarian ( 1998a, b) proposed that this fore- 
ground was electric dipole radiation from ultrasmall spin- 
ning dust grains. Although such emission from spinning 
dust had been discussed previously (see Erickson 1957, 
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Ferrara & Dettmar 1994), DL98a were the first to in- 
clude the variety of excitation and damping processes 
that are relevant for very small grains. As time went on 
alternative models for the enigmatic foreground have ap- 
peared to be inconsistent with observations!!, while the 
predictions of the spinning dust model have thus far been 
confirmed. As a result, spinning dust is now the principal 
explanation for the mysterious "Foreground X" . 

Although the model in Draine & Lazarian (1998ab) 
provided quantitative predictions consistent with obser- 
vational data, the current state of precision measurement 
of the foregrounds calls for refinement of the model, us- 
ing a better description of the complex grain dynamics 
and modifying some of the original assumptions. 

Recent studies showed that the correspondence of the 
DL98 model to observations can be improved by adjust- 
ing the parameters of the model. For instance, Doblcr 
et al. (2009) used the Wilkinson Microwave Anisotropy 
Probe (WMAP) 5 year data to show that a broad bump 
with frequency at ~ 40GHz correlated with Ha, a tracer 
of the warm ionized medium (WIM). They showed that 
this bump is consistent with predictions from a DL98 
model modified so that grains have a characteristic dipole 
moment of 3.5 D at grain size 1 nm, and the gas num- 
ber density of the WIM is jih = 0.15 cm -3 (cf. tih = 

3 For instance, the dust-correlated synchrotron emission sug- 
gested by Bennett et al. (2003) has now been ruled out (see de 
Oliveira-Costa et al. 1999, Finkbcincr, Langston, & Minter 2004, 
Boughn & Pober 2007, Gold et al. 2009, 2010). 
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0.1cm -3 in the DL98 model). 

On the theoretical front, Ali-Haimound et al. (2009) 
improved the accuracy of predictions of the model of 
emission from spinning dust using the Fokkcr-Planck 
equation for the angular velocity u). The authors quan- 
tified the deviations of the grain angular velocity distri- 
bution function from the Maxwellian approximation that 
had been used by DL98b for the sake of simplicity. With 
the other assumptions being identical to DL98b, their 
findings are not much different from DL98b's predictions. 
However, both the DL98 model and the refined model by 
Ali-Hai'moud et al. (2009) disregarded the non-sphericity 
of grains and the anisotropy in the damping and excita- 
tion processes. Obviously, this assumption is inexact for 
non-spherical grains or when there exists any anisotropy 
in the damping and excitation processes. This present 
paper is intended to go deeper into studies of grain dy- 
namics for a disk-like grain geometry, relaxing more of 
the simplifying assumptions in the original DL98b treat- 
ment. 

The main thrust of our present work is to provide a bet- 
ter description of several physical processes which have 
not been addressed in their complexity either in DL98b 
or the papers that followed. In particular: (i) the effects 
on electric dipole emission arising from the wobbling of 
the axis of major inertia of the disk-like grain around 
the angular momentum due to internal relaxation and 
the anisotropy of grain rotational damping and excita- 
tion (not yet been treated in the literature, so far as we 
are aware); and (ii) transient spin- up of very small grains 
due to single-ion collisions. 

First of all, the former process, i.e., imperfect internal 
alignment in the non-spherical grain, can alter the fre- 
quency at which the electric dipole emits. In the DL98 
model, the emission frequency is identical to the angular 
frequency lu/2tt. However, if the dipole moment is fixed 
in the grain body, then the complex motion of the grain 
axes around the angular momentum J will result in emis- 
sion at frequencies different from its angular frequency. 

In addition, imperfect alignment is essential for many 
astrophysical processes, e.g., for grain alignment (see 
Lazarian 2007 and Lazarian & Hoang 2009 for recent 
reviews). In our quest to understand the rotational dy- 
namics of grains which do not rotate about their axis of 
major inertia we capitalize on improved understanding of 
internal randomization arising from thermal fluctuations 
within the grain (Lazarian 1994, Lazarian & Roberge 
1997, Weingartner 2009, Hoang & Lazarian 2009). Dis- 
alignmcnt of the grain's principal axis from the direction 
of the angular momentum J will cause the angular veloc- 
ity to increase, leading to increased electric dipole emis- 
sion at high frequency. Our paper provides a quantitative 
description of the effect on the spinning dust emissivity. 

The latter process - collisions with ions - is impor- 
tant for small grains where the angular momentum of an 
impinging ion can be larger than the pre-collision grain 
angular momentum, resulting in a rotational excitation 
spike. In DL98b it was noted that the effect is expected 
to increase the spinning dust emissivity, but no quanti- 
tative description was given. 

While the treatment of high impulse ion collisions is 
easily performed within our Langevin code, the treat- 
ment of grain wobbling for non-spherical grains requires 
a careful and somewhat tedious modification of our treat- 



ment of the angular momentum diffusion and damping in 
the DL98 model. In particular, we have to consider sep- 
arately parallel and perpendicular contributions to grain 
damping and excitation. In addition, we provide an im- 
proved treatment of infrared emission from spinning dust 
grains. 

The structure of the paper is as follows. In §2, we 
present elements of the DL98 model and our modifica- 
tions to that model. In §3 we provide detailed calcula- 
tions for rotational damping and excitations arising from 
plasma-grain interactions and electric dipole emission for 
the disk-like grain geometry. Refined calculations for the 
effects of infrared emission are presented in §4 and 5. In 
§6, we study the effect of the grain precession on the elec- 
tric dipole emission spectrum and identify its frequency 
modes. In §7 and 8 we present our numerical techniques 
for finding the distribution functions of grain angular ve- 
locity and electric dipole emission, and benchmark cal- 
culations. In §9, we present our results for emissivitics 
from spinning dust for various idealized environments, 
and clarify the role of grain shape and differential rota- 
tional damping and excitation processes to the increase 
of peak emissivity and frequency. Discussion and sum- 
mary are presented in §10 and 11, respectively. 

2. REVISITING THE DL98 MODEL 

2.1. Elements of the DL98 model 

Here we present the grain model and our notation. 
DL98 pointed out that the abundant polycyclic aromatic 
hydrocarbon (PAH) particles required to explain the ob- 
served infrared emission provide a population of particles 
that must be spinning and emitting electric dipole radia- 
tion in rotational transitions. The smaller PAH particles 
are expected to be planar. The grain size a is defined 
as the radius of a sphere of equivalent volume. Grains 
are assumed to be disk-like with height L and radius R 
for a < ct2 and spherical for a > a.2- a 2 — 6 x 10~ 8 
cm is chosen in DL98b. The surface equivalent radius 
a s is defined to be the the radius of the sphere with the 
same surface area of the grain. The excitation equivalent 
radius a x is defined to be the radius of the sphere with 
the same J r 2 dS with r being the distance from the sur- 
face element dS to the center of mass of the grain (see 
Appendix A). 

The electric dipole moment fx of a grain arises from the 
intrinsic electric dipole moment associated with asym- 
metric molecules or substructures, and from the asym- 
metric distribution of any excess charge present. The 
latter is shown to be less important. 

A grain acquires charge through collisions with ions 
and electrons in gas, and through photocmission. As- 
suming that the charging and photocmission are in equi- 
librium (i.e., ionization equilibrium), the distribution of 
grain charge f(Z) for a given grain size can be ob- 
tained by solving the ionization equilibrium equations 
(see Draine & Sutin 1987, Weingartner & Draine 2001b). 

A grain in the gas experiences collisions with atoms 
and ions, plasma-grain interactions, infrared emission 
and electric dipole emission. All these processes result 
in damping and excitation of grain rotation. DL98b as- 
sumed that the angular velocity uj is perfectly aligned 
with the grain axis of major inertia ai , and isotropically 
oriented in space, and derived the damping and exci- 
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tation coefficients assuming perfect internal alignment. 
The dimcnsionlcss damping and excitation coefficients, 
F and G arc defined as 



F; 



r m 
3 ~ 2k B T K 



th dujj 
ujj dt 



dt 



(1) 
(2) 



where j=n, i, p and IR denote collisions of the grain 
with neutral, ion, plasma-grain interactions and infrared 
emission, (l/2)I^doJj /dt is the increase of kinetic energy 
of rotation along one axis due to the excitation process 
j; th is the damping time of the grain in a purely H I 
gas of temperature T gas , and I|| is the moment of inertia 
along the grain symmetry axis. For an uncharged grain 
in a gas of purely atomic hydrogen, Fh = Gh = I ■ 

The emissivity per H due to the electric dipole emission 
of spinning dust with angular velocity uj is 



3u 1 1 f — , dn 2 / 2mV 

— =a / da—Airuj f u 2w J - 

riH 47T tih J a da \ 3c 3 



, (3) 



where nn is the density of H nuclei, f u is the distribution 
function for the angular velocity lo, (j,± is the electric 
dipole moment perpendicular to the rotation axis, and 
dn/da is the grain size distribution function with a in 
the range from a m - m to a max . Here we take the grain 
size distribution from Draine & Li (2007), and consider 
only carbonaceous grains. In the DL98 model, for the 
sake of simplicity, f u was assumed to be a Maxwcllian 
distribution. 

2.2. Our improvements of the DL98 model 

Both the DL98 model and the refined model by Al- 
Ha'imoud et al. (2009) assumed that very small grains 
have disk-like shape and large grains are spherical, but 
they ignored the non-sphericity of grain shape when cal- 
culating the rotational distribution function and emissiv- 
ity. In the present paper, we modify the DL98 model as 
follows. 

First, we consider the rotation of a disk-like grain with 
angular velocity oj not perfectly aligned with its sym- 
metry axis (i.e. grain wobbling or imperfect alignment), 
relaxing the assumption of the perfect internal alignment 
of uj with the symmetry axis in the DL98 model. The ro- 
tational damping and excitation coefficients parallel and 
perpendicular to the grain symmetry axis resulting from 
grain infrared emission, plasma-grain interactions and 
electric dipole emission arc then derived. 

Second, we will identify the frequency modes of dipole 
emission as a result of complex motion of the grain axes 
with respect to a fixed angular momentum J. 

Third, we will find the exact distribution functions for 
angular velocity f u and electric dipole emission frequency 
f v for non-spherical grains using numerical simulations 
of the Langevin equation (LE), instead of assuming the 
Maxwellian distribution as in the DL98 model or using 
the Fokker-Planck equation (FP) as in Ali-Ha'imoud et 
al. (2009). Using the LE approach, we investigate the ef- 
fect of grain wobbling on the emission spectrum of spin- 
ning dust. The effect of internal thermal fluctuations, 
which results in deviation of the grain symmetry axis 



from the angular momentum J (Lazarian 1994; Lazarian 
& Robcrgc 1997), is also studied. 

Finally, the transient rotational spin-up due to single- 
ion collisions with very small grains is numerically stud- 
ied using the LE, in which single-ion collisions are treated 
as Poisson-distributed discrete events. 

To see how the grain wobbling can modify results from 
the DL98 model, let us consider the simple case of 3D ro- 
tation for an axisymmetric grain with principal moments 
of inertia I\ > I2 = I3 along principal axes ai, a 2 and 
a3, respectively. Denote I\\ = I\ and I± = I2 — I3. The 
ratio of moments of inertia is defined as h = 7|| / Jj_ . For 
the moment, we ignore the damping due to the electric 
dipole emission. Assuming that the rotation along three 
principal axes are independent, the mean square angular 
velocity for this case is 

(c 2 ) = (^) +2 (ci), (4) 



A) 



The mean square angular 
velocities along the parallel and perpendicular direction 
to the symmetry axis are given by (see DL98b) 



Gh 1 fc R T f 



B-tgas 



F\ 



1,-L Z ||,-L 



(5) 



where -Fji j_ and Gn j_ are total damping and excitation 
coefficients defined by equations ([1) and © correspond- 
ing to the rotation parallel and perpendicular to the grain 
symmetry axis. After some manipulations, we get 



G|| 2fc B T gas 
F\\ h 



a 
a 



1 



where 



J II F II 
Gn 



I±F± 

G± 



(6) 



(7) 



Therefore, the increase of mean square angular velocity 
becomes 



3( W f) 



)DL98 



3(w 



^,)DL98 



1 , (8) 



where (w 2 ) DL 98 = 3(w 2 ) = (G\\/F\\) X (3fc B T gas //||) . 

It can be seen that the excess of mean square angular 
velocity of the grain depends only on the "anisotropy 
ratio" 

a 11 . ( Fu G± 



r} = — - = h x 



F\G\\ 



(9) 



The anisotropy can arise from the non-sphcricity of grain 
shape (i.e, h 7^ 1) and from the differential damping and 
excitation along the direction parallel and perpendicular 
to the symmetry axis (i.e., F\\ 7^ F±, Gp 7^ G±). For 
the disk-like grain, 77 > 1 , the increase of (w ) is given 
by equation ([8]), and therefore we expect an increase of 
the frequency at the emission peak (hereafter peak fre- 
quency). Since the power radiated by a rotating grain 
is a nonlinear function of ui, we expect a substantial in- 
crease in the emissivity of spinning dust. We will study 

4 The eigenvalues of the moment of inertia tensor of the disk 
of density p, radius R, and height L are Ik = npR 4 L/2, I± = 

{(mR A L/l2)[?, + {L/R) 2 }. 
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TABLE 1 

Idealized Environments For Interstellar Matter 



Parameters 


CNM 


WNM 


WIM 


RN 


PDR 


nn (cm -3 ) 


30 


0.4 


0.1 


10 3 


10 5 


Tgas (K) 


100 


6000 


8000 


100 


1000 


T d (K) 


20 


20 


20 


10 


80 


X 


1 


1 


1 


1000 


30000 


•ni 


0.0012 


0.1 


0.99 


0.001 


0.0001 




0.0003 


0.0003 


0.001 


0.0002 


0.0002 


y = 2n(H 2 )/n H 


0. 


0. 


0. 


0.01 


0.01 



the correlation of the increase of peak frequency to the 
anisotropy in § 9.3. 

Infrequent hits by single-ion collisions are able to tran- 
siently spin-up the grain, producing spikes of angular ve- 
locity w. As a result, we expect both the peak frequency 
and total cmissivity increase when the transient spin-up 
is taken into account. This issue we address in §8. 

2.3. Idealized Environments 

Table [T] presents physical parameters for idealized 
environments where nj_ is the hydrogen number den- 
sity, T gas and Td are gas and dust temperature, x = 
Urad/wisRF is the ratio of radiation energy density u ra <j 
to the mean radiation density for the diffuse interstellar 
medium uisrf (see Mathis, Mezger, & Panagia 1983), 
n(H-2), n(H + ), n(M + ) are the molecular hydrogen den- 
sity, ion hydrogen density and ionized metal density, re- 
spectively. Physical parameters for the CNM, WNM, 
WIM, RN are similar to those in DL98b. The parame- 
ters for the PDR are taken to be similar to those inferred 
for the Orion Bar (see Allers et al. 2005). 

3. ROTATIONAL DAMPING AND EXCITATION FOR 
IMPERFECT ALIGNMENT 

Below wc present calculations of rotational damping 
and excitation for a disk-like (or cylindrical) grain with 
principal moments of inertia I» and I±. The disk has 
radius R and height L. The angular velocity u) is at an 
arbitrary angle with the grain symmetry axis ai . When 
the angular velocity uj is not aligned with ai , the rotation 
of the grain consists of the rotation about ai with angular 
velocity W|| and the rotation about an axis perpendicular 
to ai with angular velocity lu±_ . Parallel and perpendicu- 
lar components of dimensionless damping and excitation 
coefficients arc defined as in equations ((1) and @, but 
J||,u;,th are replaced by , ojii and t h n and Ij_,loj_ and 

3.1. Collisional Damping and Excitation 

For a disk-like grain, the damping times for rotation 
parallel and perpendicular to the grain symmetry axis 
due to gas collisions are given in Appendix B. For a pure 
H gas of density nn, they read 

in -1/2/30 cm -3 \ , . 

T H .|| «4.12 x 10 10 a_ 7( oT 2 1/2 ( — J T\\ s, (10) 



th._l~4.58 x 10 9 a_ 7 pT 2 1 



/2 / 30 cm' 



1 



^8 (6L\ 2/3 _ 
11 9\RJ {2L/R + 1) 



r, = - 



1/3 



r_. s, (11) 

(12) 

^) 2/3 [3+(i/i?) 2 ] xi-, (13) 



where p = p/2 gem- 3 , T 2 = T gas /100 K, a_ 7 = a/10" 7 
cm, and 
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\R 



(14) 



For dimensionless collisional damping and excitation co- 
efficients, i^n.i},!! = F {n,n,x and \\ = G {nA}:± be- 
cause they are normalized over those of purely H gas. 

It is convenient to define components of thermal angu- 
lar velocity parallel and perpendicular to the grain sym- 
metry axis at the gas temperature T gas : 



Wt. 



2k B T, 



gas 



1/2 



wt.i — 



2fc B T, 



1/2 



(15) 



3.2. Electric Dipole Damping 

The electric dipole moment fi of the grain is assumed 
to be fixed in the grain body system, and given by 



H = /iiai + /x 2 a 2 + /i 3 a 3 , 



(16) 



where fix , (12 and ^3 are components of \x along the grain 
principal axes ai , a 2 and a3 , respectively. 

The rotation of the electric dipole along ai results in 
the damping for wn, and its rotation along the axes per- 
pendicular to ai results in the damping for The 
decrease of parallel and perpendicular components of an- 
gular velocity due to electric dipole emission are given by 
(see Appendix B3) 



dt 
I±doj± 
dt 



T 2 , ,3 



1 



3fcBT gas Ted.ll 
3fcBT gas T cd ._L ' 



(17) 
(18) 



where the damping times for rotation along and perpen- 
dicular to ai are given by 



3/ 



T"ed,|| 



6^^/eBTj 



T e d,± = 



6(/i| + Ml/2)fc B T ga 



-,(19) 



and /if, = and /1 2 , = /i 2 



II — auu p 2 — p 3 — have been assumed. 

Following DL98b, the dipole moment is given by 



23 



V a 



3.8 



/3 



0.4 D 



a- 7 



Debye 2 (20) 



where (Z 2 ) is the mean square grain charge, (3 is the 
dipole moment per atom of the grain, and a_7 = 
a/10 _7 cm. Plugging equation (|2"0"|) into (|TTJ)) with the 
assumption of uniform distribution of fi along the grain 
principal axes, and using moments of inertia for the disk, 
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CNM 



wc obtain 



T ed,ll 



and 



T c d.± 



3/ 2 c 3 
Afi 2 k B T gas ' 

1.6 x 10 n p 2 a- 7 x 
3/ 2 c 3 



(i?/L) 4 / 3 



[(a x /a) 2 (Z 2 )+3.8(/?/0.4D) 2 a_ 7H2 



:4.0 x 10 w p 2 a 8 _ 7 x 



(fi/L) 4 / 3 (l + l(£/i?) 2 ) 2 
[(a x /a) 2 (Z 2 ) +3.8(/3/0.4D) 2 a_ 



For grains larger than 122, T od || = T c d,_i_. 

3.3. Damping and Excitation by the Plasma 

The problem of plasma-grain interactions for perfect 
internal alignment was studied in DL98b, and refined 
by Ali-Hahnoud et al. (2009). Here we consider the 
damping and excitation by the plasma for the case of 
imperfect internal alignment. We assume for simplicity 
that /j, is directed along a 2 axis. We consider only neutral 
grain (Z g = 0, as in DL98b), and assume a spherical 
grain with a cylindrical excitation equivalent radius a cx 
(e.g. Ali-Hahnoud et al. 2009). 

For an ionized gas with ion density nj, ion mass mi, 
and charge Z- u the dimensionless excitation coefficients 
parallel and perpendicular to the grain symmetry axis 
from plasma drag are given by 




G n II I — — 



"H 



TO; 



1/2 



\k B T ( 



cx 



du g 



l,-L> 



(23) 



different vaiues of cj||/aj Ti || and ui±/ojt,± m the CNM. Electric 
dipole moment /3q = 0.4 D is adopted. 

in the ultraviolet. At wavelengths 3000 < A < 1500 A, 
photoabsorption is primarily due to 7r — > n* electronic 
transitions, excited by electric fields parallel to the plane 
defined by the C atoms. PAHs - especially PAH ions - 
are also able to absorb at longer wavelengths, although 
the absorption is much weaker. The electronic transi- 
tions responsible for the absorption are not certain, but 
we will here assume that this absorption also involves 
transitions with electric dipole matrix elements parallel 
to the C atom plane. 

Suppose the interstellar radiation field to consist of a 
unidirectional component providing a fraction 7 of the 
local energy density of starlight, plus an isotropic com- 
ponent. Let n* be a unit vector parallel to the direction 
of propagation of the unidirectional component of the 
starlight. Let the dimensionless factor 



where e is the elementary charge, u = v/vt, with v being 



the velocity and 


14 = 


2fc B T' gas /mi, 
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/a cx ' V 
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\2k B T gas 



uQu 
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1/2 



(24) 



where I\\ x (uf2|i_i_//) are given by equations (|C12j) and 
(|C13j) (see Appendix C2 for the definition of b max and 
detailed calculations). Using the Fluctuation-Dissipation 
theorem (see Lazarian 1995; Lazarian & Roberge 1997 for 
the application to grain diffusion coefficients) , we obtain 

Figure [1] presents G p .|| and C7 Pj j_ for different values 
W||/ojt,|| ancl w±/ w t.± for CNM. It can be seen that G Pi || 
and G p _± increase as the grain rotates slower, and for 
wiij_ < 10 -3 , they converge to the value for non-rotating 



gram, i.e, w^j. 



0. 



4. PHOTOEXCITATION OF PAHS 

Small PAHs arc planar molecules, consisting of sp 2 - 
bonded C atoms in a two-dimensional hexagonal struc- 
ture, with peripheric H atoms. PAHs absorb strongly 



U 



"MMP 



(25) 



where u* is the starlight energy density, and mmmp = 
6.85 x 10~ 14 ergs cm -3 is the estimate of Mathis, Mezger 
& Panagia (1983, hereafter MMP) for the local starlight 
energy density. Let N a bs,o be the photon absorption rate 
for the PAH if illuminated by isotropic starlight with U = 
1. Assuming the unidirectional and isotropic components 
to have similar spectra, the PAH photoexcitation rate is 
(Sironi & Draine 2009) 

N abs = UV(6)N abSt0 , 



(26) 

*(f?)= (l- 7 )+rj(l + cos 2 fl) , (27) 

cosfJ^n* • ai, (28) 
where ai is a unit vector normal to the C atom plane. 

5. REVISED ROTATIONAL EXCITATION AND DAMPING 
BY IR EMISSION 

DL98b discuss the various processes contributing to 
the rotational excitation and dcexcitation of very small 
dust grain. The dominant processes include direct colli- 
sions of atoms and ions with the dust grain, torques on 
the dust grain due to electric fields produced by passing 
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ions ( "plasma drag" ) , and exchange of angular momen- 
tum with the electromagnetic field via absorption and 
emission of photons. Because a dust grain radiates many 
more photons than it absorbs, the emission processes are 
much more important than the change of angular mo- 
mentum resulting from photon absorptions. 

If the grain is rotating, then there will be a tendency 
for the infrared emission to, on average, reduce the grain 
angular momentum ( Martin 1972). This process was in- 
cluded in the study by DL98b. As noted by Ali-Hai'moud 
et al. (2009), the DL98b expression for the rotational 
damping torque from infrared emission was too small by 
a factor of two owing to an algebraic error. Here we 
make a more refined estimate of rates for excitation and 
deexcitation of PAH rotation due to infrared emission. 
We employ the PAH model of Draine & Li (2007, here- 
after DL07); this grain model reproduces various obser- 
vations of infrared emission, and therefore should provide 
an improved estimate of rotational excitation and damp- 
ing from infrared emission. The present treatment also 
takes into account anisotropy in the infrared emission 
from planar PAHs. 

Table 2 lists the character assumed for each of the vi- 
brational modes in the DL07 PAH model. 



TABLE 2 

Optically- Active Vibrational Modes 



X(fim) 


in-plane fraction 


identification* 


3.30 


1 


C-H stretch 


5.27 





oop C-H bend overtone? 


5.70 





oop C-H bend overtone? 


6.22 


1 


C-C stretch 


6.69 


1 


C-C stretch 


7.417 


1 


C-C stretch 


7.598 


1 


C-C stretch 


7.850 


1 


C-C stretch 


8.33 


1 


ip C-H bend 


8.61 


1 


ip C-H bend 


10.68 





oop mono C-H bend 


11.23 





oop mono C-H bend 


11.33 





oop mono C-H bend 


11.99 





oop duo C-H bend 


12.62 





oop trio C-H bend 


12.69 





oop trio C-H bend 


13.48 





oop quartet C-H bend 


14.19 





oop quartet C-H bend 


15.90 


2/3 


? 


16.45 


2/3 


? 


17.04 


2/3 


? 


17.375 


2/3 


? 


17.87 


2/3 


? 


18.92 


2/3 


? 


15. 


2/3 


FIR continuum? 


* ip = 


in-plane mode; oop 


= out-of-plane mode 



Having in mind PAHs, we consider a grain with dy- 
namical symmetry, with eigenvalues I\ > 1% = I3 of the 
moment of inertia tensor. Let ai,a2 and a3 be unit vec- 
tors parallel to the principal axes corresponding to mo- 
ments of inertia I\ , I2 and I3 , respectively. The vibra- 
tional modes will be assumed to have oscillations that 
are either parallel to ai ( "out-of-plane" ) or perpendicu- 
lar to ai ( "in-plane" ) . Let the angular momentum vector 
have components J|| parallel to ai, and J± perpendicular 




A(yu.m) 

Fig. 2. — Moments of the infrared emission from a PAH ion with 
effective radius a = 6.31 A (Nq = 105 carbon atoms), heated by 
the average interstellar radiation field (U = 1). The radiated power 
(top panel) and photon emission rate (middle panel) are dominated 
by 5 ^; A ;$ 15/^m emssion, but the rotational damping coefficient 
(bottom panel) has an appreciable contribution from emission out 
to A 200/jm. The sawtooth discontinuities in the spectra at 
A H 50/xm arise from approximations in the calculational method 
(see text). 



to ai (i.e., in the a2-a3 plane). 

Suppose that the grain is absorbing photons at a rate 
N a b s , and radiating photons with a power per unit fre- 
quency P v . For U < 10 4 , small PAHs are in the single- 
photon-heating regime, with virtually all of the energy 
of an absorbed starlight photon being radiated before 
the next photon absorption occurs. This allows us to 
assume a universal time-averaged emission spectrum for 
each PAH, with the emission intensity proportional to 
the photon absorption rate: 

P m ,u « U^{9)P m , u , , (29) 

where -P, rl; „,o is the time-averaged power per unit fre- 
quency radiated by in-plane (to = ip) or out-of-plane 
(m = oop) vibrations of a single PAH molecule or ion 
heated by isotropic radiation with U = 1. Figure [2] 
shows vP mM fi for a PAH with Ac = 105 C atoms, where 
Pm.v,o is obtained by solving the equations of statisti- 
cal equilibrium for a PAH subject to stochastic heat- 
ing (Draine & Li 2001, Li & Draine 2001b), using the 
"thermal-discrete" method, where the emission from a 
PAH at a given energy is estimated using a thermal ap- 
proximation, but the downward transitions are treated as 
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I0 1 



Fig. 3. — Rate of photon emission TV (upper panel) and (A) (lower 
panel) as a function of Nq, the number of C atoms in the PAH, 
for PAHs in the local interstellar radiation field. 
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IR Rotational Damping times 

A 




10 2 



•t ± (PAH + 
t ± (PAH°) 



10 3 



10 4 



N„ 



Fig. 4. — Time scale for damping of angular momentum parallel 
to the symmetry axis (Jy) and perpendicular to the symmetry axis 

(Jx)- 



discrete transitions to a lower energy level. The model 
employs a heat capacity based on a realistic spectrum 
of vibrational modes (Draine & Li 2001), and updated 
vibrational band strengths (Draine & Li 2007) which, 
for a suitable PAH size distribution, closely reproduce 
observed IR spectra. The sawtooth discontinuities in 
the model spectrum arises from approximations in the 
treatment. The long- wavelength emission from a real 
PAH would presumably be dominated by a finite number 
of emission lines, but the continuous spectrum obtained 
here is expected to provide a reasonable representation 
for the distribution of emitted power over wavelength. 
The photon emission rate 



where 



N„ 



-Pm,v.odv. 



(30) 



(31) 



The middle panel of Figure [5] shows the integrand in eq. 
(|3"Tj) for the time-averaged rate of photon emission N m 
by in-plane and out-of-planc modes, for PAH neutrals 
and ions. 

The rotational damping also depends on the photon- 
weighted wavelength, 



1 



A 



(32) 



for m = ip or m = oop. The lower panel of Figure 
[5] shows the integrand in eq. (f3"2")l . For the PAH emis- 
sion properties assumed by Draine & Li (2007), the small 
fraction of the total power radiated at A ^ lOO^rn con- 
tributes substantially to the integral in eq. (|32|) . but does 
not dominate the integral. 



Figure [3] shows how N m and (A) depend on the the 
number Nq of C atoms in the PAH. 

Let J be the grain angular momentum. We assume 
that the grain has no intrinsic helicity, so that a nonrotat- 
ing grain is equally likely to emit left- or right-circularly 
polarized photons. If the grain is initially non-rotating, 
IR emission will result ir0 



=5* 



N, 



oop,0 



N, 



bs,0 



N, 



ip.O 



(33) 
2fr 2 .(34) 



The infrared emission will, on average, remove angular 
momentum from the grain, with rates (sec Appendix B2) 



dt 
dJ ± 
~dt 



in 



in 



3^ 
27r/||C 



t/*(#)A lp ,o(A) 4p J|| 
1 



2tt/i c 



C 



(35) 
.7436) 



The characteristic time scales for rotational damping 
by IR emission are 



-J\\ 



1 2tt7|C 



7"IR,|| 
7"IR,X 



(dJ l{ /dt) m ~ UV{0) 3h 

-J ± 1 2tt7 ± c 

(dJ±/dt) m ~ U^{6) 3H 

x (l/2)iV ip , (A) 4p + A oop , (A) oop 



,(37) 



(38) 



The rotational damping times tjru and t\r i ± are 

5 The factor 2 in front of h 2 is pointed out by Silsbce et al. 
(2010). 
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shown in Figure 0] for PAHs in the average interstellar ra- 
diation field. We show separate damping times for PAH 
and PAH + , but the timescale for changes in PAH ioniza- 
tion are short compared to the time scale for changes in 
angular momentum, so one should take an appropriate 
weighted mean: 

T m,\\ = Mrji ( PAH °) + (1 - /oh^u (PAH+), (39) 
Tm,± -M R a( PAH °) + (1 " /o)r IR ^(PAH+), (40) 

where /o is the neutral fraction for PAHs of that size. 

Also shown is the rotational damping time estimated 
by DL98b for graphitic grains. Our new estimates for 
Tirii and T[r.j_ are smaller than the DL98b estimate by 
factors of 1-6, depending on size, charge state, and the 
orientation of ai relative to the angular momentum J. A 
factor of 2 is attributable to an algebraic error in DL98b 
(see Ali-Hai'moud ct al. 2009); the remaining differences 
are due to use of realistic PAH emission properties rather 
than the power-law heat capacity and power-law opacity 
used by DL98b, as well as different assumptions regard- 
ing the grain moment of inertia. 

Using equations (TTJ) and ([2]), the parallel and perpen- 
dicular components of the dimensionless coefficients of 
damping and excitation from infrared emission are then 
given by 



Gm ,|| 



T H,|| _ 

tir,|| ' 
dJ 2 



-Fir, j_ — 



T~m,± 



(41) 



dt '2/||fc B T gas 
t;U^{0) ' r ,, „ r H ,||, 



Gm,± — 



2I\\kBT gas 

r H,_L 



(42) 



,dJ 2 



dt 2I±k B T gas 



N oop .o + hiNabs.O + ^Vip,o) 



2h 2 



2I±kBT gas 



-r H ,x(43) 



where 7|| and Jj_ are given by equation 

The new coefficients for rotational excitation and 
damping by infrared emission derived in this section are 
compared with those from the DL98b model in Figures [S] 
and [6] for WIM, RN and PDR. There the electric dipole 
damping rate for parallel and perpendicular direction is 
also shown. For spherical grains (i.e. a > 6 x 10 _8 cm), 
the parallel and perpendicular components of -Fir and 
Gir are not the same as the electric dipole damping rate, 
i.e., there exists the anisotropy in the damping and ex- 
citation. Such a high anisotropy of the damping and ex- 
citation from infrared emission can result in important 
effects of the dynamics of spinning dust grains. 



6. 



EFFECT OF GRAIN PRECESSION ON THE DIPOLE 
EMISSION 



Earlier studies (DL98ab, Ali-Haimoud et al. 2009) as- 
sumed that the frequency of dipole emission from rotat- 
ing grain is the same as the rotational frequency u/2tt. 
Since the electric dipole is fixed in the body system, the 
rotation and precession of the grain with respect to fixed 
J can induce some modification for the dipole emission 
frequency from its rotational frequency. 
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Fig. 5. — Comparison of our new damping and excitation coeffi- 
cients arising from infrared emission, -Fir, Gir with those in DL98b 
for WIM with /o ~ 0. The ratio of damping times th/t b< i for ro- 
tation parallel and perpendicular to the grain symmetry axis also 
shown. 
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Fig. 6.— Similar to figure [5] but for RN and PDR. Here the 
neutral fraction of PAH /o = 1. 

Let us consider the simple case of a disk-like grain 
where the dipole lies in the plane &i&2'- 



(44) 



in the grain body system. The torque-free motion of 
this grain with angular momentum J and ratio of inertia 
moments h consists of the precession of the symmetry 
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Fig. 7. — Torque-free motion of the grain in an external system 
with z|| J and x, y _L J described by three Euler angles 8, <j> and 

axis ai about J with constant angle 9 and rate <j>, and 
the rotation of grain itself about the symmetry axis with 
the rate ip where <f>, ip and 9 are Euler angles (see Figure 
[7]). The precession and rotation rates are respectively 
given by (Landau & Lifshitz 1976) 



— , ip = (1-h)— cost 



(45) 



Precession and rotation of the grain with respect to J 
results in an acceleration for the dipole moment: 



H = /i||ai + ^j_a 2 , 



(46) 



where ai and a 2 arc functions of time and given in Ap- 
pendix F. 

The instantaneous emission power of the dipole mo- 
ment is defined by 



P ed (J,9,t) 



2 

3<? 



Ml 



(47) 



The emission power of the grain is obtained by aver- 
aging (|47|) over the torque- free motion of the grain, i.e., 
over angles ip and <p in the range from to 2tt. The final 
result for the case = and fi±_ = \fj,\ is 



Ped(J,0) _ 1 + COS 2 (9 

2\i\ /3c 3 ~ 2 



(<t> + ^Y 



i ip + 0-0 3 )(l - cos I 



) 2 + sin 2 i 



'•(48) 



To find the frequency of dipole emission, we perform 
the Fourier transform for the components of dipole accel- 
eration fl. Normalized square amplitude of the Fourier 
transform for the components jl x (or jl y ) and jl z are 
shown in Figure [8] for 9 = 15 and 40° for a disk-like 
grain with h = 1.5. We can see that the emission spec- 
trum from the electric dipole corresponds to three fre- 
quency modes, ip/2ir, resulting from the component jl z , 
and (<p±ip)/(2n), arising from the component jl x (prjly). 



The emission power from the former mode is negligible 
compared to that from the later modes, which have fre- 
quencies given by 



>±i> J hT (h- 1)| cos I 



2tt 



'II 



2 71 



(49) 



in which the major mode (<j> + ip)/2ir for 9 < 90° (or 
(4> — f/')/27r for 9 > 90°) is much stronger than the sec- 
ond mode (4>—tP)/2tt ((<p+7p)/2Tr). When 9 increases, the 
amplitude of higher frequency mode increases, and that 
mode becomes important as 9 approaches 90° (see Fig. 
[5]). However, for the fast internal relaxation, the angle 9 
fluctuates with small amplitude about J, and the lower 
frequency mode is dominant. Therefore, we are inter- 
ested only in the emission of this lower frequency mode 
in the present paper. The ratio of emission frequency 
of the dominant mode to the rotation frequency w/27t is 
given by 



4>- |V>| _ h- (h- 1)| cos 0| 



w/2tt 



h 2 i 



< 1. 



(50) 



It turns out that the dominant emission frequency is 
smaller than the rotation frequency, and it is very close 
to the rotation frequency for very small 9. The cmissivity 
at the frequency v per H is given by 



Ju_ 

n H 



1 1 

47T 7lH 



dn 
da—jva 
da 



(51) 



where 



j, >a = p e d(J,e)f v (J,6), 



(52) 

and f v (J,9)dv is the probability of the dipole emission 
in [v, v + dv\ . 

If we assume that all of the emission is at the domi- 
nant frequency given by ([50]) , then equation (|52|) can be 
rewritten in terms of integrals over J and 9 by introduc- 
ing a Delta function: 



d9dJP ed {J,9)S(v' -v)fjf g . 



(53) 



where v' = v'(J, 9) is the emission frequency correspond- 
ing to J and 9, fjdJ is the probability of the grain having 
angular momentum in [J, J + dJ], and fgdO is the prob- 
ability of the angle being in [9, 9 + d6]. 

It is possible to calculate j v>a using (|53|) . However, 
we can assume that the emission power is approximately 
given by a power law of the emission frequency 



2 2 
3c33 



Mi(2^) 4 



(54) 



We found that the emissivity obtained using (|5 1 1) and 
(|52[) with this approximation does not differ more than 
10% from more rigorous calculations using (fpU]). There- 
fore, in the following, we calculate the emissivity using 
equations (|5H52[) with the approximated emission power 



P e d{v) given by equation (fM]) . 

7. NUMERICAL METHODS 

The most important element in calculating the emis- 
sivity from spinning dust grains is to find the probability 
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Fig. 8. — Normalized emission spectrum for a disk-like grain 
(h = 1.5) with /j. 1 1 = 0, for two values of the angle 9 between 
the principal axis ai and the angular momentum J. The two com- 
ponents of \FT(il x )\ 2 (or |FT(/i y )| 2 ) arc indicated by circles, while 
the component of |FT(/t z )| 2 is indicated by a triangle. 

distribution f u and f v for rotation resulting in emission 
at frequency v for a given grain size a and a given electric 
dipolc moment fx. Here we present a numerical approach 
using the Langevin equation to obtain those distribution 
functions. 

7.1. The Langevin Equation 

The Langevin equation is a well-known equation for 
describing stochastic dynamics (see Gardner 1983). In 
terms of its assumptions and applicability it is equiva- 
lent to the well known Fokker-Planck equation. In par- 
ticular, the equation can handle only situations when the 
increments of the variable are smaller than the variable. 
Thus, the Langevin equation cannot be applied directly 
for treating the high impulse collisions of spinning dust 
grains with ions. Therefore in what follows we use a hy- 
brid approach, describing the grain dynamics for most 
of the computational time using the Langevin equation 
and describing the infrequent (but high impact) collisions 
with ions as impulsive events. 

The angular momenta of the impacting ions are drawn 
from distribution functions that include focusing by 
the attractive potential when approaching neutral or 
negatively-charged grains. For describing grain dynam- 
ics the Langevin equation was first suggested in Robcrge, 
DcGraff & Flaherty (1993) and found more applications 



in later works on grain alignment (see Roberge & Lazar- 
ian 1999, Hoang & Lazarian 2008, 2009). 

7.1.1. ID rotation 

In an idealized situation when the grain rotates only 
about its axis of major inertia ai, namely ID rotation, 
the rotation of the grain is determined by the value of 
angular velocity uj = W||, and distribution function /„ = 

2tt/w0 Note that for this case, tjj = Thji , and cjt = ^t,|| i 
and notations || and _L are omitted in this section. 

For a given grain size, the evolution of w can be de- 
scribed by a stochastic differential (Langevin) equation. 
In dimensionless units of lu' = uj/ljt and t' = t/ra, the 
Langevin equation reads 



duj' 



F tnt co' 



2w M r H 



3r c d 



dt' 



-dq, (55) 



where dt' = dt/rn is the time step, dq is a random vari- 
able with variance (dq 2 ) = dt, and the total damping 
Ftot and diffusion coefficients Gtot are given by 



Ftot — F\ + F n + Fir + F p , 
Gtot = G\ + G n + Gir, + G p , 



(56) 
(57) 



where Fj and Gj with j=i, n, IR, p arc damping and ex- 
citation coefficients corresponding to ion-grain, neutral- 
grain interactions, infrared emission and plasma drag. 

To solve equation (|55p for a/, we numerically integrate 
it for N time steps with total time T. After each time 
step dt, the value u/ is updated via equation 



J i+1 



FtoM 



2^ 3 T H 



3r, 



cd 



dt' 



G, 



-dq, (58) 



for i running from 1 to N. The initial value lu' q is a 
random variable in the range [0,1] generated from a uni- 
form distribution function^ The obtained solution is 

binned and used to find the distribution function fu . 
Wc note that fu is f° r the rotation along one axis only. 
By assuming isotropic distribution of the vector u>, the 

distribution function f u oc u; 2 /^ 1 ' with normalization 
condition L fudui = 1. Results for fu are presented 
in §7.3. 

7.1.2. Grain Wobbling: 3D rotation 

When the angular momentum is not completely 
aligned with the grain axis of major inertia (i.e., 3D ro- 
tation), its orientation within the grain body has to be 
accounted for in finding f u and /„. 

Problems of internal alignment, i.e. to what extent the 
angular momentum J and the grain axis of major inertia 
are aligned have been always at the focus of the quantita- 
tive treatment of grain alignment (see Lazarian 2007, for 
a review). The initial works (e.g. Jones & Spitzer 1967) 
usually assumed a Maxwcllian distribution for J. Then 
Purcell (1979) noticed that the internal relaxation of en- 
ergy within a wobbling interstellar grain may be fast (i.e, 

6 This situation is similar that in the DL98 model and in the 
improved treatment by Ali-Hai'moud et al. (2009). 

7 The initial value is unimportant as it is "forgotten" if the 
simulation is sufficiently long. 
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the timescale of internal relaxation is smaller compared 
to the gas clamping time j3- Due to the fast internal relax- 
ation demonstrated in Purcell (1979) for a period of time 
the works on grain dynamics assumed that the grain axis 
of major inertia must always be directed along the grain 
angular momentum, as this alignment minimizes the en- 
ergy of the rotating grain. However, Lazarian (1994) 
noticed that this cannot be true for thermally rotating 
grains and restored the notion of grain wobbling as an 
essential element of grain dynamics. 

In the 3D rotation case, we consider the evolution of J 
in the lab coordinate system instead of u>. For a general 
grain shape, J is completely determined by three compo- 
nents J x , y ,z in the inertia coordinate system eie 2 e 3 (see 
Appendix E). The evolution of J in time is also described 
by three Langevin equations 

dJj = Ajdt + \J Badqi, for i = x, y, z, (59) 

where dqi is the random Gaussian variables with (dqf) = 
dt, At = (AJi/At) and B vl = ((A J ?; ) 2 /At) are damp- 
ing and diffusion coefficients defined in the inertial co- 
ordinate system. In dimensionless units J' = J/Jii^t ii, 
t' = t/ra n, equation (|59[) becomes 



'2\ 



where (dq i 
A' = -J'K 



dt', and 



(60) 



tot, | 
7' 3 



7 H sm 



(cos 4 + 7 edsin 4 0- 

6 T od 



h 3 + 3h 



Bu 



2In knT, 



- T H,|| 



shVflcos^ (61) 



(62) 



where 



' cd 



r cd,|| 
T H,|| 



7h 



-Ftot,J_TH, 
-Ptot,|| T H.J 



7ed = 



I|| r cd,|| 
I±T e d,l 



,(63) 



and 9 is the angle between ai and J, -Ftot,|| j -Ftot,_L are 
total damping coefficients parallel and perpendicular to 
the grain symmetry axis, and Bu are given in Appendix 
E. 

7.2. Effect of Internal Thermal Fluctuations 

Very little is known about internal relaxation in mi- 
croscopic spinning grains. Thus, in what follows, for our 
description of grain wobbling we consider two extreme 
models. The first model assumes fast internal relaxation, 
with the distribution of deviations of the axis of grain 
maximal inertia from the direction of J determined by 
the dimensionless ratio J 2 / 'Ink^T^, where Td is the tem- 
perature of the grain body (Lazarian & Roberge 1997). 
The second model assumes no internal relaxation, with 
a Maxwellian distribution function for the angle between 
ai with J. 

8 Purcell (1979) introduced a new magneto-mechanical effect, 
which he termed Barnett relaxation, which he showed to be able 
to align 0.1 fim grains over a time scale of the order of a year. Later 
another magneto- mechanical effect, termed nuclear relaxation, was 
introduced in Lazarian & Draine (1999). Nuclear relaxation is 
much faster than Barnett relaxation for 0.1 (im grains, but is not 
applicable to the tiny spinning grains considered here. 



For a given value of angular momentum, and in the 
presence of fast internal relaxation, the angle 9 between 
the grain axis of major inertia ai and J follows the local 
thermal equilibrium distribution function 



/lte(0) = A exp 



J 



2/|ifcB7d 



(1 + [h - l]sin 2 



■ A exp — J 



'2 J gas 



(1 + [h - l]sin 2 0) , (64) 



where h = Iu/Ix, J' = J/Ihuj^u . and A is the normaliza- 
tion constant such that J* j 'lte(#)27t sin 9d9 = 1 (Lazar- 
ian & Roberge 1997). 

When the timescale of internal relaxation is much 
longer than that of gas damping time, i.e., without inter- 
nal relaxation, the angle distribution function is simply 
Maxwellian: 



fuw(9) 



h 



4tt 



hs'm 



,3/2 



(65) 



where /m w (#)27r sin 6>d6> = 1 (see Jones & Spitzcr 
1967; Lazarian & Roberge 1997). 

To follow the evolution of J in time in the presence 
of internal relaxation, we use the following algorithm. 
First, we generate initial values at t = for J', angle £ 
between J and ei and azimuthal angle x from uniform 
distribution functions. Thus, J' z (t = 0) = J' cos£, J' x {t = 
0) = J' sin £ cos \ and J' y (t = 0) = J' sin £ sin \- Then we 
solve equations (|60|) with the time step dt to get J' x , J' y 
and J' z and then J' at the moment t+dt. Due to the inter- 
nal relaxation, the angle 9 at that moment is unknown. 
Therefore, at the end of each timestep, 9 is assumed to 
be a random angle generated from the distribution func- 
tion (|64[ fast internal relaxation) or (|65l without internal 
relaxation). With 9 known, we can obtain the value of 

angular velocity &i = ^/o; 2 + ui'j^ with its components 

W|| = (J//||)cos0, and uj± = (J/ 1±) sm9. The dipolc 
emission frequency of the dominant mode Vi is calcu- 
lated by (|49|) . This process is repeated for N timcsteps. 
The obtained value of angular velocity and dipole emis- 
sion frequency at the end of each time step Ui and i>i are 
binned and will be used to find the distribution function 
f u and /„. 

7.3. Benchmark Calculations 

To benchmark our Langevin code, we first run sim- 
ulations for the perfect internal alignment case to find 
f u = (l/27r) /„ to compare with the analytical distri- 
bution function obtained from the Fokker-Planck (FP) 
equation (Ali-Hai'moud et al. 2009). Physical parame- 
ters for different phases are given in Table 1. For very 
small PAHs, the dust temperature fluctuates due to tran- 
sient heating by single photon. Therefore, we assume a 
smaller dust temperature T d = 10K for the CNM, WNM 
and WIM, and T d = 20 and 40K for the RN and PDR, 
respectively, for smallest grains with a < 7 x 10 _8 cm. 
Rotational damping and excitation coefficients in Fig- 
ure (j2"0"]) are adopted for the benchmarks. We note that 
the results of LE simulations depend on the choice of 
right time step. Very small timesteps are inefficient, 
while large timesteps can introduce error. We define the 
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Fig. 9. — Comparison of normalized uj 3 fui obtained from analyt- 
ical solution of the FP equation and from numerical simulations of 
the Langevin equation for grain sizes a = 3.56 and 4.5 A. 
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Fig. 10. — Comparison of peak frequency with and without elec- 
tric dipole damping (ed) obtained from simulations of LE and FP 
for WIM (upper panel) and PDR (lower panel). 



timestep by dt 
show that e = 10" 



e min[l/F tot , l/G to t, T H /r od ]. Our tests 
1 is a good choice. We run the LE 
simulations for N = 10 7 steps. 

Figure [9] shows normalized w 3 / u for a grain with elec- 
tric dipole moment parameter /3 = fa and sizes a = 3.56 
and 4.5 A obtained from the FP equation (solid lines) 
and the LE simulations (dashed lines) for the WIM0 
It can be seen that the LE simulations are in excellent 
agreement with analytical solutions of the FP equation. 

Another important parameter we need to benchmark 
is the peak frequency as a function of grain size. For the 

9 Hereafter, the distribution function f u for the case /? = 0q is 
presented. 



perfect alignment case, the peak frequency is given by 
equation (160) in Ali-Hai'moud et al. (2009). Using the 
distribution f u from LE simulations, we can obtain the 
peak frequency by maximizing the function uj f u . Re- 
sults from these two approaches are presented in Figure 
[TU1 for the WIM and PDR for both cases with and with- 
out electric dipole damping. In all cases, results from 
LE simulations coincide with those from the FP. For the 
PDR, the results with and without electric dipole damp- 
ing are not much different because the electric dipole 
damping is negligible. 

8. IMPULSIVE EXCITATION BY SINGLE-ION COLLISIONS 

DL98b showed that for grains smaller than 7 x 10 _8 cm, 
the angular impulse due to an individual ion-grain col- 
lision may be comparable to the angular momentum of 
the grain. Thus, infrequent hits of ions can result in 
transient rotational excitation for small grains. 



Let 



be the mean rate of ion collisions with the 



grain, given by 

-/(Z g = 0)n i7 ra 2 



'icoll ' 



/(Z g )n i7 ra 2 

z e ^o 



m;7T 



1/2 



1/2 



-<I> 



+ 



Z g Zie 2 



,(66) 



where $ = (2Z 2 e 2 /ak B T gas ) 1/2 , g(x) = 1 - x for x < 
and g(x) = e~ x for x > 0, and f(Z s ) is the grain charge 
distribution function (see Appendix D). The probability 
of the next collision occurring in [t, t + dt] is 



dP 



i dt. (67) 

The rms angular momentum per ion collision (S J 2 ) is in- 
ferred by dividing the total rms angular momentum to 
the collision rate, and its final formula is given in Ap- 
pendix D. 

The upper panel in Figure [11] presents the collision 
rate, Tj~ oll compared to the rate of electric dipole damping 
for the grain with /3 = 0.4D, r^ 1 , for the CNM, WNM 
and WIM. The effect of single-ion collisions is determined 
by the critical size a cr i corresponding to t^ o11 = . 
For the WIM, it can be seen that for grains with a < 
a cr i = 8.6 x 10 _8 cm, the electric dipole damping time is 
shorter than the time between two ion collisions. The 
active range of single-ion collisions for the WNM is a < 
a cr i = 1.2 x 10~ 7 cm. The critical size is smaller for the 
CNM with a cri = 3.8 x 10~ 8 cm. Although, the single-ion 
collisions are important for the CNM, WNM and WIM, 
they are not important for the RN and PDR. 

To account for single-ion collisions, we first use the 
Poisson distribution to generate the time intervals be- 
tween successive collisions ATicoiij, which is a random 
variable in this case for a given Ti co ii. Thus, the nth col- 
lision occurs at the moment ticoii.n = A'Hcoiij. We 
then run the integration for the LE (i.e., eq. l60| over 
time to find J and u. The nth collision deposits an an- 
gular momentum AJ n that is assumed to be randomly- 
oriented, with magnitude drawn from the distribution 
function with the rms {5 J 2 ) appropriate for the grain 
charge and radius. The impinging ions are taken to be 
protons for the WNM and WIM, and a mix of protons 
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Fig. 11. — Upper panel: Comparison of ion collision rate Tj n 
with the electric dipole damping rate T~ d of grains with electric 
dipole moment parameter /3 = 0.4D for the CNM, WNM and WIM. 
The filled circle symbols denote the critical size a cr ; when T.Zi, = 
r^j . Lower panel: The evolution of the angular velocity w/tJ T ii 
as a function of t/ T H ,|| f° r the case with and without impulses from 
single-ion collisions in the WIM for a = 4A. 

and C + ions, with fractional abundances xh and xm from 
Table 1, for the CNM. For the WIM and WNM in which 
the focusing effect due to negatively-charged grain is neg- 
ligible, we can assume that the grain charge is given by 
the mean charge (Z). But for the CNM in which the 
focusing effect is important, we consider the grain in var- 
ious charge states and find the corresponding f u and 

fr- 

To identify the role of single-ion collision on the distri- 
bution function, we first consider the perfect alignment 
case and solving ID LE (Equation (f5"5j0 . The resulting 
ia for the WIM is shown in the lower panel of Figure [IT] 
for a grain size a = 4A. It can be seen that due to the 
electric dipole damping, the grain rotates subthcrmally 
most of the time with u < Slo. Ionic impulses result in 
the transient rotational spin-up followed by the continu- 
ous damping by electric dipole emission. 

Distribution functions obtained from LE simulations 
in the presence of impulses and the grain wobbling are 
compared with those when only the grain wobbling is 
considered (i.e., without impulses) in Figure [12] for the 
WNM (upper panel) and the WIM (lower panel) and for 
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Fig. 12. — Normalized w 4 /^ in the presence of impulses from 
single-ion collisions for the WNM (upper panel) and the WIM 
(lower panel) for two typical grain sizes. Impulses broaden the 
emission spectra to higher frequency (see diamond symbols). 

two grain sizes 3.55 and 4.0A. It can be seen that the 
ionic impulses extend the distribution function to higher 
frequency, but change only slightly the peak frequency of 
the spectrum. This effect is more important for smaller 
grains. 

9. EMISSIVITY 

Here we calculate the emissivity at the frequency v 
using equation (|51|) with f u obtained from the LE sim- 
ulations and the emission power P e< x(y) from equation 
(|54]l. To compare with earlier studies, we also calculate 
the emissivity using equation ([3]) with the Maxwellian 
distribution function f u from DL98b and the analytical 
distribution function from the FP equation (FP). 

The emissivity is calculated assuming that 25% of the 
grains have the electric dipole moment parameter /3 = 
2/3 , 50% have /3 = P and 25% have /3 = 0.5/3 with 
Po = 0.4D. The electric dipole moment fi is assumed 
to be isotropically distributed along the grain principal 
axes. We adopt the models of grain size distribution from 
Draine & Li (2007) with the total to selective extinction 
i?v = 3.1 and the total carbon abundance per hydrogen 
nucleus be = 5.5 x 10 -5 for diffuse environments CNM, 
WNM and WIM, and i? v = 5.5, b c = 2.8 x 10~ 5 for the 
RN and PDR for carbonaceous grains with a m ; n = 3.55 
A and a max = 100 A. 

9.1. Grain Wobbling: Imperfect Internal Alignment 
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Here we present our new calculations of emissivity ac- 
counting for imperfect internal alignment for the ideal- 
ized environments listed in Table 1. We run simulations 
of LE (cq. [60]) to find the distribution functions /„ for 128 
grain sizes in the range from a m ; n to a max . For plasma 
drag and infrared emission, we use our new damping and 
excitation coefficients from Figures ((T|) and ([5]). For col- 
lisional damping and excitation, we use the results in 
DL98b as shown in Figure (j20|) . 



9.1.1. Distribution function f u 

For a more instructive comparison of our results with 
earlier studies, here we present the distribution func- 
tion of rotational frequency f u instead of the distribu- 
tion function of dipolc emission frequency /„, which is 
slightly different from the former one for the grains of 
the planar geometry. 

Figure [T3l shows normalized w 4 / w for the case of grain 
wobbling with and without fast internal relaxation, com- 
pared to that obtained using the FP method in Ali- 
Hai'moud et al. (2009) for a grain size a = 4 A for the 
WIM and PDR. To find the distribution function for the 
FP case, we use the parallel damping and excitation co- 
efficients for the distribution function given by equation 
(34) in Ali-Hai'moud et al. (2009) because the grain is 
assumed to rotate along the symmetry axis. 

It can be seen that in both environments, the grain 
wobbling results in a higher peak frequency, and the en- 
tire distribution functions are shifted to higher frequen- 
cies. Furthermore, the distribution function for the case 
without internal relaxation is extended to higher fre- 
quency than that with fast internal relaxation because 
the higher rotation frequency in the later case corre- 
sponds to higher internal dissipation rate of energy (see 
dotted line in the lower panel). 

9.1.2. Emissivity 

We calculate the emissivity per H using the distribu- 
tion function /„ found in the previous subsection for dif- 
ferent environments in the presence of fast internal re- 
laxation and without internal relaxation. 

In Figure Q3J we compare our results in the presence of 
fast internal relaxation with those obtained using meth- 
ods in DL98b and Ali-Hai'moud et al. (2009) for perfect 
internal alignment for the WIM and RN. The parallel 
damping and excitation coefficients are taken in use to 
find the distribution function for the later cases. 

First, it can be seen that the grain wobbling results 
in increases of the peak frequency. Compared to re- 
sults from the FP equation method in Ali-Hafhoud et al. 
(2009), we see that the peak frequency increases from 
~ 23 GHz to - 33GHz for the WIM and from - 106 to 
~ 161 GHz for PDR. Also, the peak emissivity of spin- 
ning dust in the later case is increased by factors of ~ 1.5 
for the WIM and ~ 1.8 for PDR, respectively. 

The emissivity per H for various idealized environ- 
ments are summarized in Figure[15]with and without fast 
internal relaxation, respectively. The solid lines denote 
the emission spectra for the case of fast internal relax- 
ation. The shaded areas represent the expected spectra 
when the internal relaxation varies from fast to very slow. 
The subplot represents j v / jv,FV with j v and j„.fp being 
the emissivity at the peak of the emission spectrum for 
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Fig. 13. — Comparison of normalized uj^fui for the case of perfect 
internal alignment using the FP method (dashed line, FP) in Ali- 
Hai'moud et al. (2009) with that in the presence for grain wobbling 
obtained from simulations of LE, both with fast internal relaxation 
(LE, F) and without internal relaxation (LE, N). Shading areas 
denote the transition from fast internal relaxation (solid line) to 
without internal relaxation (dot line). For both the WIM and 
PDR, the peak frequency of distribution function increases, and 
the spectrum shifts to higher frequency in the later cases. 
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Fig. 14. — Comparison of emissivity spectrum for the grain wob- 
bling (solid line) with those from the ID rotation case (dot and 
dashed lines) and the fast internal relaxation (LE, F) is taken 
into account. The peak frequency of WIM increases from ~ 23 
to ~ 33GHz and its increases from ~ 106 to ~ 161 GHz for PDR, 
relative to that from the FP equation method in Ali-Hai'moud et 
al. (2009). Emissivity increased by a factor of ~ 1.5 and 1.8 for 
WIM and PDR, respectively. 
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Fig. 15. — Emissivity per H in the case of grain wobbling for 
various environments. Gray shaded areas represent the transition 
between fast internal relaxation (F, solid lines) to no internal relax- 
ation (N, dot lines). Subplot shows jv/j v ,FP with j^ ; fp being the 
peak emissivity from the method of FP equation, as a function of 
the increase of peak frequency Af pea k for fast internal relaxation 
(F) and no internal relaxation (N). 

the imperfect alignment case and that from the FP equa- 
tion method, as a function of the increase in the peak 
frequency = ^ pca k — ^pcak.FP of emission spectrum. 

Figure [T5] shows that the increase of emissivity and 
peak frequency for the grain wobbling case is important. 
The peak frequency is increased from 8 to 55GHz (i.e., by 
factors of 1.2 to 1.8) different environments (see subplot 
in Figure I15[) . The peak emissivity is increased by fac- 
tors 1.47 (1.42), 1.5 (1.4), 1.52 (1.37), 1.8 (1.2) for WNM, 
WIM, CNM and PDR, respectively, corresponding to the 
cases with fast internal relaxation (without internal re- 
laxation), except for the RN where the peak emissivity 
enhancement can reach a factor of 4.1(2.2) (see the sub- 
plot). In addition, the enhancement of peak frequency 
for the case without internal relaxation is smaller than 
that for fast internal relaxation. That can be explained 
by the fact that some emissivity from its peak can be 
transferred to higher frequency, resulting more extended 
spectra present in the case without internal relaxation 
(dot line in Figure ITSl) . 

The highest increase of peak emissivity present in 
the RN seems difficult to explain. We note that for 
this medium, the infrared damping and excitation is 
the most important process. Therefore, in addition to 
the anisotropy due to the non-spherical geometry, the 
anisotropy in the damping and excitation due to IR emis- 
sion results in such a substantial increase of peak emis- 
sivity and peak frequency (see equation [8]) . The corre- 
lation of the increase of peak frequency and emissivity 
with the anisotropy for the grain wobbling case will be 
investigated in §9.3. 

9.2. Emissivity in the Presence of Single-Ion Collisions 

To see the effect of single-ion collision on the emissivity 
of spinning dust, consider first the case of perfect align- 
ment. We are interested only in the CNM, WNM and 
WIM where the single-ion collisions are important. 

We run the LE simulations for the grain wobbling tak- 
ing into account the impulses for grain size from a min to 
a cr ; for these environments. We consider both the fast 
internal relaxation and without internal relaxation. The 



resultant emissivity is compared with that from the case 
without impulses in Figure [16] for the WNM and WIM. 

It can be seen that the ionic impulses broaden effi- 
ciently the emission spectra to the higher frequency, but 
they change slightly the peak frequency. For the WIM, 
the net effect of impulses and grain wobbling enhances 
the peak emissivity by a factor of 1.73, and the peak 
frequency is increased from ~ 23 to ~ 35 GHz. The 
corresponding increase in peak emissivity for the WNM 
is a factor of 1.58, and the peak frequency is increased 
from ~ 19 to ~ 27 GHz. By subtracting the effect of 
grain wobbling, we see that the impulses can increase 
the peak emissivity by ~ 23% and 11% for the WIM and 
WNM, respectively, but they increase slightly the peak 
frequency (see Figures [15] and [T6|) . 

For the CNM, the lower panel of Figure [T7l shows that 
the effect of impulses on the total emissivity is weaker 
than for the WIM and WNM. The reasons for that are as 
follows. First, although the focusing effect is important 
for the grain in the negatively charged state, resulting in 
rotational impulses at high frequency (see the dot line in 
the upper panel of Figure[T7|). the probability of the grain 
in this negative charge state is relatively small. Second, 
since the effective range of impulses of the CNM is much 
narrower than that for the WIM and WNM (see the up- 
per panel of Figure []TJ , its emissivity, that is obtained 
by integrating over the entire grain size distribution, is 
obviously less affected by impulses compared to the WIM 
and WNM. 

We note that for a small grain in the CNM, WNM and 
WIM, during a short interval of several ti co ii, the grain 
experiences some rotational spikes due to single-ion colli- 
sions as shown in the left lower subplot of Figure [TT] The 
peak frequency and width of the spectrum are related to 
the mean square value, (w 2 ), which is averaged over total 
integration time T. Therefore, the overall effect of the 
transient spin-up is to broaden the emission spectrum 
and increase the total emissivity as seen in Figure \W[ 

9.3. Effect of the non- sphericity and anisotropy in 
rotational damping and excitation 

To understand why the grain wobbling induces the sub- 
stantial increase of peak frequency and emissivity found 
in the previous sections, let us consider a disk-like grain 
shape and when internal relaxation is important. 

To investigate the effect of grain wobbling for both 
spherical grains (a > 02) and disk-like grain (a < 02), 
we run LE simulations to find f u and f pea k for the case 
of imperfect internal alignment assuming fast internal 
relaxation, for grain size from a m - ln to a max 

We consider here the RN and PDR for which electric 
dipole damping is subdominant. Results are shown in 
Figure [15] where solid and dashed lines represent i/ pca k 
obtained using the DL98 model and LE simulations, re- 
spectively. The ratio of f pca k to that of DL98 model is 
depicted in the upper right subplot ™ The excitation 
anisotropy ratio a\\/a± (see eq. [9]) is shown in the lower 
left subplot. 

First of all, for the PDR, the peak frequency from LE 
simulations is slightly higher than that from the DL98 
model for spherical grains (a > a-2 = 6 x 10~ 8 cm). For 

10 When the electric dipole damping is negligible, the distribu- 
tion f u is Maxwellian, therefore i/ pca k,FP = fp C ak,DL98 
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Fig. 16. — Emissivity per H obtained for the WNM (upper panel) 
and WIM (lower panel) without ionic impulses using the FP, and 
with impulses using our LE simulations for grain wobbling. The 
spectra arc efficiently broadened as a result of impulses. Fast in- 
ternal relaxation is more efficient at highest frequency, resulting in 
narrower spectra compared to those for the case without internal 
relaxation (red color line). 

very small disk-like grains (a < 6 x 10 _8 cm), the peak 
frequency is increased by a mean factor of ~ 1.5. For the 
RN, one interesting feature is that the peak frequency 
from LE is much larger that that from the DL98 model 
for spherical grains(cf. to the PDR). We attribute the in- 
crease to the anisotropy in the grain rotational damping 
and excitation due to the infrared emission, which is the 
most important process of grain excitation for the RN. 
For spherical grains in the range from a2 to 10~ 7 cm, 77 
increase from 1.2 to 1.4 for the PDR, and from 1.8 to 2.3 
for the RN (see subplots in the lower left corner of Figure 
fl~8|) . Hence, the higher anisotropy in the RN corresponds 
to the higher increase of peak frequency compared to the 
PDR case. 

To see clearly the dependence of the increase in i/ pca k 
in the grain wobbling case, we plot t'peak/^'peak.DLgs as a 
function of 77 for the PDR in Figure HH The solid line 
corresponds to 



(68) 




^pcak,DL98 



where diamonds show f pea k from the LE simulations. We 
see that the ratio t'peak/^peak.DLgs increases with 77 as 
expected. For small 77, ^ poa k from the LE simulations is 
similar to that from equation (|68l) . and some difference 
appears when r\ increases (see Figure [T8|) . This may stem 
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Fig. 17. — Effect of impulses due to single-ion collisions for the 
CNM: normalized ui 4 /,(, with f' w is the normalized distribution 
function, as a function of ui 1 = u/ojt,|| where Z is the grain charge 
Z (upper panel) and emissivity (lower panel). Gray and blue 
shaded areas correspond to the transition between fast internal re- 
laxation (solid line) and without internal relaxation (dot line) for 
the case without and with impulses, respectively. The focusing ef- 
fect for the grain with charge Z = — 1 enhances the high-frequency 
tail (dot line in upper panel) 

from two reasons. First, the angular velocity distribution 
function for ?y ~ 1 is Maxwellian, but it differs from 
the Maxwellian distribution when 77 increases. Second, 
the assumption that u> x ,y,z are independent in deriving 
equation (|68p is not valid because they have complex 
motion in the inertial coordinate system. 

In other words, the anisotropy due to non-sphericity of 
grain shapes and the anisotropy in the grain rotational 
damping and excitation both act to increase the peak 
frequency of the emission. 

Finally, we note that the total emissivity due to spin- 
ning dust scales as w 4 / u , thus, the increase of ^ pC ak nat- 
urally is accompanied by the increase of the peak emis- 
sivity as observed in Figure [TBI 

10. DISCUSSION 

10.1. Comparison with earlier studies 

The DL98 model was the first proposal to explain the 
20 — 40 GHz anomalous foreground emission in terms of 
rotational emission from ultra-small dust particles. Nu- 
merous observations (Dobler et al. 2009 and references 
therein) have confirmed the DL98 model predictions. 
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Fig. 18. — Peak frequency as a functions of grain size for the PDR 
(a) and RN (b). Dashed and solid lines denote the result obtained 
from LE simulations ^ pca k] and from the DL98 model ^peak,DL98: 
respectively. The ratio fpeak/^peak.DLQS versus a is shown in the 
upper corner subplot, and the anisotropy r\ = a^/a± is shown in 
the left lower corner. 
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Fig. 19. — The ratio of peak frequency t / pe ak/'''peak,DL98 as a 
function of a nisotropy r\ obtained from LE compared to the ideal 
case (eq. I68j l for the PDR. The gap between 77 = 1.4 and 1.8 arises 
from the change in the grain geometry from sphere to disklike. 



In the DL98 model the calculations were mostly an- 
alytical. The present study used a different approach, 
namely, the Langevin equation. We successfully bench- 
marked our calculations against the Fokkcr-Plank calcu- 
lation in Ali-Hai'moud et al. (2009) for the perfect align- 
ment case, and then introduced additional effects not 
included in DL98. In particular, our approach allows 
us to treat discrete high impact impulses arising from 
single-ion collisions with the grain. Most importantly, it 
enables us to study the spinning dust emission for the 
case of imperfect internal alignment with and without 
fast internal relaxation. 

Among grain rotational damping and excitation pro- 
cesses, plasma drag is shown to be important for Molec- 
ular Clouds and CNM. Ragot (2002) considered the in- 
teractions of the plasma with spinning dust from the 
point of view of generation of waves. We believe that 
for the case we consider, the dust-plasma interaction is 
accurately approximated by interaction of the grain with 
individual, uncorrelated, passing ions. 

DL98b noted that the typical rotational quantum 
number of even a small PAH would be J/h^i 10 2 , al- 
lowing the dynamics to be treated classically. Ysard 
& Verstraete (2009) formulated the problem quantum- 
mcchanically; the close correspondence between their re- 
sults and those of DL98b confirms that a classical treat- 
ment is entirely adequate 03 

While our studies confirm the general validity of the 
DL98 model of spinning dust emission, they also show 
ways of making the model predictions more accurate. 
For instance, our refined treatment of grain dynamics 
provides peak emissivity predictions which differ by fac- 
tors from ~ 2 to 4 from the original ones, depending on 
the environments. These differences may be detectable 
with high precision future observations. 

The largest uncertainty within the model is the value 
of the dipole electric moment. The distribution of dipolc 
moments adopted in DL98b was only an educated guess. 
Moreover, this distribution could be affected by inter- 
stellar processes and change from one media to another. 
For instance, Dobler ct al. (2009) showed the anomalous 
emission in the 5-year WMAP data exhibits a bump at 
~ 40 GHz. This peak frequency is larger than the predic- 
tion by the original DL98b model but could be explained 
if a modified distribution of electric dipole moments is 
adopted (Dobler et al. 2009), but the new effects consid- 
ered here - misaligned rotation and impulsive excitation 
by ions may also explain the observed emission spectrum. 

10.2. Effect of single-ion collisions 

DL98b pointed out that for the CNM, WNM and WIM 
the angular momentum transferred in a single ion colli- 
sion can be larger than the angular momentum of the 
rotating grain, but neither DL98b or the subsequent pa- 
per by Ali-Hai'moud et al. (2009) treated the impulsive 
nature of the ion collisions. 

In the present paper, we describe ion-grain collisions 
as Poisson-distributed discrete events. When an ion hits 
the grain surface, the grain undergoes transient spin-up 
and then gets damped, mainly by electric dipole damp- 

11 In discussing some effects, e.g., dust-plasma interactions, 
DL98b did include quantum limitations for the angular momen- 
tum transfer, i.e. that A J cannot be less than h. 
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ing. Due to strong electric dipole damping, small grains 
spend most of their time rotating subthcrmally between 
two rotational spikes due to single-ion collisions. Ionic 
impulses are shown to be important for grains smaller 
that - 8.6 x 10" 8 , 1.2 x 1(T 7 and 3.8 x lCT 8 cm for 
the WIM, WNM and CNM, respectively. Our quan- 
titative simulations show that the impulses extend the 
high-frequency tail of the emission spectrum but change 
the peak frequency only slightly. They can enhance the 
peak emissivity by 23% and 11% for the WIM and WNM, 
respectively. 



10.3. Effect of Grain Wobbling 

The DL98b model assumed perfect alignment of the 
grain axis of maximal moment of inertia with the in- 
stantaneous angular momentum. This assumption is 
valid only for strong internal relaxation and suprather- 
mal (much faster than thermal) rotation of grains (Lazar- 
ian 1994). The tiny grains we deal with rotate ther- 
mally or subthcrmally. For such a rotation, the grain 
axes are not coupled with the angular momentum. Due 
to imperfect internal alignment (i.e., grain wobbling), 
the anisotropy arising from grain shape, and differential 
damping and excitation (e.g. due to infrared emission 
and plasma drag) obviously affect the grain mean rota- 
tion rate and therefore the emissivity. Increased grain 
anisotropy, as measured by the anisotropy ratio r/ (eq. [9]) 
produces an increase in the peak frequency and emissiv- 
ity. 

As we do not know the rates for internal relaxation 
in microscopic molecular size grains, we considered two 
cases: very rapid internal relaxation, and very slow inter- 
nal relaxation. Using the Langevin equation, we showed 
that for the former case, the peak emissivity of spinning 
dust increases by a factor of ~ 2, and the peak frequency 
of the spectrum increases by factors from 1.4 to 1.8, de- 
pending on the environment. Grain wobbling has a larger 
effect in broadening the emission spectrum when the in- 
ternal relaxation is slow. Especially, for the RN when the 
anisotropy of damping and excitation from infrared emis- 
sion is very high, the peak emissivity can be increased by 
a factor of 4 and peak frequency is increased by a factor 
of 2. 

For the WIM, the joint effect of grain wobbling and 
impulses due to single-ion collisions causes the peak fre- 
quency to shift from ~ 23 to ~ 35 GHz. It is inter- 
esting to note that the new predicted peak frequency is 
close to the bump with the anomalous emission found in 
the WMAP five year data (Dobler et al. 2009). In their 
paper, using the DL98b model, they explained the ob- 
served spectrum with a modified model of spinning dust 
in which the grains have broader distribution of electric 
dipole moments. They also varied the number density of 
the WIM to produce spectra with higher peak frequency. 
These modifications do not seem to be necessary with our 
improved treatment of grain dynamics. 



11. SUMMARY 

The present paper provides a refinement of the DL98 
model of spinning dust emission. Our principal results 
are outlined below: 

1. We derived the rotational damping and excitation 
coefficients arising from plasma drag, infrared emission, 
gas drag and the electric dipole damping for small disk- 
like grains with imperfect internal alignment. 

2. A nonspherical grain precessing around its angular 
momentum J radiates primarily at two frequency modes. 
The lower frequency is dominant and is not the same as 
the rotational frequency oj/2it. 

3. Using the Langevin equation, we calculated the 
emissivity of grains with disaligncd rotation, both with 
and without fast internal relaxation. For fast internal re- 
laxation, allowing for disaligncd rotation, the peak emis- 
sivity of spinning dust is increased by a factor ~ 2 for 
the CNM, WNM, WIM and PDR; and by a factor of - 
4 for the RN, compared to those from the FP equation 
method. Fast internal relaxation also shifts the emis- 
sion spectra to higher frequencies factors of 1.4 to 2 for 
these environments. Slow internal relaxation broadens 
the emission spectrum and reduces the enhancement in 
the peak emissivity. 

4. The increase of peak emissivity and frequency 
of spinning dust in the grain wobbling case compared 
to those from the DL98 model has two sources: the 
anisotropy of the moment of inertia tensor, and the 
anisotropy in the rotational damping and excitation pro- 
cesses (e.g. infrared emission and plasma drag). Higher 
anisotropy results in higher enhancement of peak fre- 
quency and emissivity. 

5. We devised a method to include the high- impact ion 
collisions within our statistical treatment of grain rota- 
tion. We found that for the CNM, WNM and WIM, 
single-ion impulses result in strong broadening of the 
emission spectrum to higher frequency and increases the 
emissivity of spinning dust. Impulsive excitation is most 
efficient for the WIM, but it is not important for the RN 
and PDR. In all idealized environments, the transient 
spin-up is subdominant compared to grain wobbling in 
terms of enhancement of emissivity, but it is dominant 
in spectrum broadening for the CNM, WNM and WIM. 

6. For the WIM, the net effect of grain wobbling and 
transient spin-up from single-ion collision increases the 
peak emissivity by a factor ~ 1.7, and the peak frequency 
increases from ~ 23 to ~ 35 GHz. That increase in peak 
frequency brings the model into agreement with the 5- 
year WMAP data without modifying the electric dipole 
momentum distribution or number density of the WIM 
as was proposed by Dobler ct al. (2009). 
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APPENDIX 
A. GRAIN PROPERTIES 
Al. Grain shape and size 

Here we present grain model and assumptions, notations we use throughout the paper. Grains are assumed to be 
disklike with radius R and height L for size a < and spherical for a > a^- We chose — 6 x 1CP 8 cm as in DL98b. 
The "surface-equivalent" radius a s and the "excitation equivalent" radius a x of the grain were defined in DL98b. 

Ali-Ha'imound et al. (2009) defined the radius of cylinder with the same J r 2 sin 9 dS with rsin6* is the distance from 
the grain surface to the symmetry axis, thus 

4ttc4 = \j> r 2 &m 2 6dS. (Al) 

The expressions for a s , a x and a cx are given in DL98b and Ali-Hai'moud et al. (2009). We adopt the models of grain size 
distribution from Draine & Li (2007) with the total to selective extinction Ry = 3.1 and the total carbon abundance 
per hydrogen nucleus be = 5.5 x 10 -5 for diffuse environments CNM, WNM and WIM, and i?v = 5.5, 6c = 2.8 x 10~ 5 
for the RN and PDR for carbonaceous grains with a m ; n = 3.55 A and a max = 100 A. 

A2. Moments of Inertia 

Consider a PAH with Ac carbon atoms, plus peripheral H atoms that make a minimal contribution to the total mass. 
We take the surface density of C atoms in a planar PAH to be that in a single plane of graphite: oc = 3.7 x 10 15 C cm -2 . 
We will assume that PAHs with Ac < N\ consist of a monolayer, but larger PAHs will be multilayered, with number 
density of C atoms nc = 1-1 X 10 23 cm~ 3 . We approximate the PAH as a disk of radius R and height L, with 



1 + A max(0, A* /3 - A, 173 )] , R 2 = (A2) 



nc 

Neglecting the contribution of the H atoms, the moments of inertia are 

R 2 (R 2 L 2 \ 

7||=7Vcm c — , I± = A c m c f — + — j , (A3) 

Wc take Ni = 10 2 and A = 0.4; with this choice, the PAH consists of two layers when Ac = 364. 

B. ROTATIONAL DAMPING 

Gas-grain interactions, plasma-grain interactions, infrared emission, and radio emission damp the rotation of small 
grains. In this section we present the rates for these processes. 

We consider a disk-like grain with temperature Td in gas of temperature T gas , H nucleon density n-n = n(H) + 2n(H2), 
and with n(Hc) = O.lnjj. Following DL98, we "normalize" the various drag processes to the drag which would be 
produced by "sticky" collisions in a pure H gas of density njj: thus for drag process j, wc define the dimcnsionless 
quantity FjA\,± corresponding to the rotational damping parallel and perpendicular to the symmetry axis ai, such that 
the contribution of process j to the drag torque is 

For the disk- like grain in a pure hydrogen gas of density nn, Fj.\\.± = lj and the damping times for rotation along 
and perpendicular to its symmetry axis are given by 

3/|| 3-Ti. 
TH '" 47ra4 x n H mH(2fcBT gas /7rm H ) 1/2 ' TH-± 7ri? 4 nHTOH(2fc B Tg as /7rm H ) 1/2 cj_L ' 

where 

1 /L\ 3 L 1 f L\ 2 1 . x 

^=6UJ + R + 2\r) + 2- (B3) 

Using standard parameters for the interstellar diffuse medium (ISM) and moments of inertia 7ii , 7j_ for the disk, we 
get 



TH,[| ~ 4 ' 12 x 10 10 pa_ 7 T- 1/2 ( 30 ° m 3 ) r,| s, and r H ,T w 4.58 x 10 9 pa_ 7 T 2 " 1/2 ( 3 ° Cm - ) T ± s, (B4) 



where p = p/2 gcm- 3 ,T 2 = Tg as /100 K, a_ 7 = a/10" 7 cm, where a = (3/4) 1 / 3 L 1 / 3 J R 2 / 3 , and 

8 {6L\ 2/3 1 /4\ 1/3 /i\ 2/3 r , T/m2l 1 



t » = 9{-r) wirtty r - = UJ [r) [ 3+ ^ x ii (B5) 
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We will use Th,||,_L as a fiducial time scales for the different sources of rotational damping. Thus the rotational 
damping times due to process j are 

(B6) 



The linear drag processes are additive: 



^11, -L = ^n,ll,-L + ^i,||,-L + ^P,|| 



IR.II 



(B7) 



where F n , F[, F p , and Fir are the contributions from neutral impacts, ion impacts, plasma drag, and thermal emission 
of infrared photons. 

B.l. Plasma Drag 

For a grain with an electric dipole moment /x, the interaction of this electric dipole moment with passing ions in the 
plasma results in a damping and excitation (Anderson & Watson 1993). The damping due to this process is obtained 
using the Fluctuation Dissipation theorem, i.e., F p = Cp. || , J_ with Gpj| t j_ are derived in Appendix C2. 

B.2. Rotational Damping by IR Emission 
Let a 2 , a 3 , ai be axes fixed in the grain. Consider an oscillator with dipole moment 

/! = /i 2 a 2 coscjgi = [(a 2 cosw i + a 3 sinwoi) + (&2 cosw i — a 3 sinwoi)] • 

Now consider rotation of the grain around ai with angular velocity uj r : 

a 2 = x cos co r t + y sin ui r t, a 3 = y cos ui r t — x sin w r t, 
where x, y, z define an inertial coordinate system. The dipole moment is 

H = — [xcos(w r + u>o)t + y sin(w r + Wo)i + xcos(w r — uj )t + y sin(w r — u>a)t] , 



(B8) 
(B9) 
(BIO) 



i.e., a dipole (// 2 /2) rotating with frequency (wo+w r ), plus a dipole (/z 2 /2) rotating retrograde with frequency (u)Q—0J r ). 
The net loss of angular momentum is 



2( M2 /2) 2 ( Ur + W ) 4 2( M2 /2) 2 (w - LO r ) 4 



3c 2 



3c 2 



LJQ — 0J r 



4^ 2 w 2 



(Bll) 



The same result for the torque is obtained for rotation around a 3 , but rotation of this oscillator around a 2 would not 
result in a torque. 
Extending this to three incoherent oscillators: 



H = cos w 02 i + /i 3 a 3 coswo 3 £ + //iai coswoii 



wc find 



dt 



(B12) 



(B13) 



where i, j and k run from 1 to 3 so that if i = 2 then j = 3 and k = 1 and so on; and replace w i ~ w 02 ~ W03 by Wo- 
lf we now consider a grain with symmetry axis ai and /U2 = A*3- The power radiated by in-plane and out-of-planc 
vibrations is 



E 



ip 



NiphujQ = + ^3)3^31 E oop = N oop hwu = [i\ ^-3, 



where Ni p and N oop are the rates of photon emission by in-plane and out-of-plane vibrations. Then 



d 
dt 
d_ 
dl 



J± 



\2ncI\\J 



N ip \) J,, 



N IP X + N oop X J ± 



(B14) 

(B15) 
(B16) 



where J± = J 2 , A = 2irc/u>o is the emitted wavelength. 
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B3. Electric Dipole Damping 

The spinning grain will radiate power, which results in rotational damping of grain. The associated rotational 
damping time depends on u>. 
In the grain body system, the dipole is given by 

H = /iiai + /i 2 a 2 + ^3&3- ( B1 7) 

For simplicity, we assume that /i 2 pa ss = ■ The angular velocity us in general is not aligned with ai- Thus, 
the rotation of the grain consists of the rotation about ai with uiu and the rotation with lo± about axes perpendicular 
to ai. The former rotation of the dipole results in the damping for wm, and the later results in the damping for co±. 
Consider first the rotation with angular velocity wii about ai. After a time t, the dipole moment becomes 

fi = /i 2 (cosu;||ta 2 + sinaj||ia 3 ) + /x 3 (— sinu;||ta 2 + cosa;||ia3), 

= fi2V2 — a2sin (J^ +v\\tj + ascos (J^ + co\\tj , (B18) 

The increase of angular momentum of the grain resulting from the rotation is then 

dt ' 3c 3L ^"~ ^ J ~ 9c 3 ' 

where we assumed = i/x 2 . Hence, 



dJ 2 , , An 1 



j. dw\\ 4/z 2 



dt ■ 9c3-f- ^ B2 °) 



Similarly, for the rotation with lo±_ about a 2 and a 3 axes, we get 

dui ± 4/i 2 3 



Since a> is not parallel to ai, let us define characteristic damping times associated with and uj±, corresponding 
to the rotation along and perpendicular to ai- 
As in DL98, we can define the damping times due to electric dipole emissions as 



1 duJ \\,±\ %-L W j|,_L 1 

'||,± dt J ed 3fc B T gas r cd || ± 



(B22) 



Then, combining equations (|B20j) , (|B21I) and (|B22|) . we obtain 

3/ 2 c 3 (R/L) 4 / 3 

T ^W = a 2i. hr = L6 x 1Ql V 2 a?.7 x r , / y2\ \ L L /n A VF s ' ( B23 ) 

4/x 2 fc B T gas [(a x /a) 2 (Z 2 ) + 3.8(/3/0.4 D) 2 a_ 7 JT 2 

and 

^ - VfceT^ - 4 '° X 10 ' a " 7 X [(a x /a) 2 <Z 2 ) + 3.8(0/0.4 D) 2 a_ 7 ]T 2 S ' (B24) 

where ^ 2 from DL98b (eq. |2"0")) has been used and (Z 2 ) is the mean square grain charge. 

Using the distribution charge function f(Z) and mean charge (Z) obtained by solving equation of ionization equations 
(see DL98b), we compute damping coefficients F n , i 7 !, F p , F\r and r e d- Our new results for F p , Fir and r e d are presented 
in Figures [TJ [5] and El respectively. For comparison, in Figures [20] we present results from the DL98 model for the 
WIM and RN. 

C. ROTATIONAL EXCITATION 
An initially stationary grain will have its rotational kinetic energy increasing at an average rate 

j t Qj||,±«f, ± ) = x (G ntl± + C^^ + G p ,|| )± + G ir ,||,_l), (CI) 

where the normalized excitation rate G n is due to impacting neutrals, Gj is due to impacting ions, G p is due to plasma 
drag, and Gir, is due to recoil from infrared emission, and ||, _L denote excitation along and perpendicular to the grain 
symmetry axis. 

CI. Recoil from Thermal Collisions and Evaporation 
For collisional excitations, similar to collisional drag, we have G n » = G n ,i_ and G,- n = Gj,j_. 
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Fig. 20. — Left: Damping coefficients F n , Fx, Fp, Fm an d T e c j/tn and excitation coefficients G n ,Gj,G p and Gir from DL98b model as 
functions of grain size a for WIM. Right: Similar to Figure l20l but for the RN. 



C2. Excitation by the Plasma 
Consider a disk-like grain with dipole moment fi, assumed to be parallel to a.2: 

**=Haa. (C2) 

Let us define the inertial coordinate system x, y, z. Grain is assumed to be rotating with angular velocity wy and 
parallel and perpendicular to the symmetry axis. Initially, a2||x, a3||y and ai||z. At the time t, the orientation of the 
grain in the lab system can be determined by two Euler angles 9 and (f> where 9 and 4> are the rotation angles of the 
grain about ai and a3, respectively^ They are given by 

ai = cos 9 z + sin 9 cos </> x + sin 9 sin <fi y, (C3) 
a2 = — sin 9 z + cos 9 cos <fi x + cos 9 sin <fi y, (C4) 
a3 = cos cf> y — sin <fi x. (C5) 

The trajectory plane of the incoming ion of charge Z\e with e being the electron charge is yz where z is parallel to 
the initial direction of ion. The electric field produced by the ion at the grain position is 

E = — = ^-(cosa y + sma z). (Co) 
The torque induced by the dipolc-clcctric field interaction is given by 

— = [fi x E] = (x[sint? cos a + cos (9 sin0 sin a]) (C7) 

+ (— cos 6* cos0 sinay + cos^ cos</> sinaz), 
where equation (|B17|) has been used. By projecting (|C7|) onto two rotational axes ai and a3, we obtain 
du\\ -Zie\n\ duj ± -Z { e\^\ . . . 

— — = — — cos0 cos ot, — — = — — (cost/ sin a + sin sm<p cos a), (Co) 

da in ov da l±bv 

where bv = r 2 da/dt has been taken. 

The mean rotational excitation due to one ion is then given by 

«-> 2 > - (W 1 - «^ ,3> " (^) V > (C9) 

where the factor of 2 comes from the fact that a in the range from to tt/2, instead from — ir/2 to tt/2, or time t from 

12 Ali-Hai'moud ct al. (2009) assumes that initial grain orien- a j. r\ ■ ■ 1 -tu a j. n u 4 -4-1, 4- 

, . . t>o,0o- Our assumption is special with On = <t>n = 0, but without 

tation is described by On, tpn, and later, they average over random , ■ , . , a ■ , 

J vrui 1 j loosing physical enccts 
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to oo, and 

ft \ 2 ft \ 2 

(CIO) 



Z|| = dacosift cosa^ , X± = (^J da(cosd sina + sinfl sine 



For a neutral grain and the trajectory of ion is straight line, the polar angle a and t are related to each other as 
follows: 

vt 6 1 vdt 

sing = =. cos q = — =, da = ; — tttoT— r-i (Cll) 

Vb 2 + v 2 t 2 Vb 2 + v 2 t 2 [l + (^/6) 2 ] b y ' 

where we assume that t = corresponding to the time when the ion is closest to the grain, i.e., r = b. We have, at the 
time t=0, the direction of the grain is random, then at t, we have <f> — 4>q + ui\\t, 9 = 9 a + uj±t with </> , 9q are random 
variables. 

Therefore. Xu and X\ become 



= x z K 1 (x) z , (C12) 



, n ( f i ,\ff , \ f f°° , cos(xr) 

Xii (a:) = / cosffl cos a da) = / coswiii cos a da = / ar — — 

" U / U " / Uo (l + r 2 ) 3 /2 

X±(y,y±) = (^j [cosw^t sina + sinw^t sinaj||t cos a] da^j , 

= T^tS + - V-Ki(y-)j = U(y) + 1 -[ V+ K l {y + ) - y.K^y.)^ , (C13) 

where x = uj\\b/v,y = Lo±b/v and y± = (u>u ± ujj_)b/v. K n (x) is the modified Bessel function of the second kind, and 
H(y) is given by 

To calculate (|C14|) numerically , we note that when the grain rotates much faster than the ion motion, i.e., y = 
Lv±b/v >> 1, the torque gets averaged out. And a truncation of (|C14[) is y = 1.3 is adopted without introducing large 
uncertainty. 

The total excitation rate from all ions with density n\ following normal distribution of velocity is given by 



In equations (|C15[) and (|C16[) . the value 6 max for a neutral grain is given by 



'max — ^cx 



<t> 2Z 2 e 2 v 
Ocx\/l + -, $ 2 = — > « = —with v T = (2fc B T gas /m i ) 1 / 2 , (C17) 



^cx^ri- 1 gas 

and for a charged grain, 6 max takes the form (see Spitzer 1941) 

Z,Z-,e 





f 1 _2^4) 1/2 for2¥ i > -- 2 



a n x 1 



and the upper limit of the integration over the impact factor is constrained by the Debye length Ad = 

(kT gas /(4nn e e 2 )) 1/2 . 
The dimensionlcss excitation coefficients for grain-plasma interactions arc then 

d((^||) 2 ) rf((<W) 2 ) JxT g , ± 

Gp <»- eft 2fc B T gas ' Gp < X ~ <K 2fc B T gas ' [LjU) 
Plugging (IC15I) and (|C16|) into the above equations, we obtain 

/ \ 1/2 / ry \ 2 /o 



G *=^UJ Ura U."- (C20 > 
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where 

Using the Fluctuation-Dissipation theorem, we obtain F p ^ = G Pi ||, and F Pt ± = G Pi j_. In Figure [TJ we show our new 
results for G p n and G Pi ± as functions of grain size for neutral grains in the CNM. 

Excitation coefficients from DL98b G n , G;, G p and Gir are also shown in the lower panels of Figure l20l for the WIM 
and RN, respectively. 

C3. Rotational Excitation by IR Photons 

Consider the case of a grain with dynamical symmetry, with principal values of the moment of inertia tensor 
I\ > I2 = I3, where I\ is the moment of inertia for rotations around the ai axis. We assume the emission to come 
from optically active vibrations that are either "in-plane" or "out-of-plane" , and we assume the "in plane" oscillators 
to be symmetrically distributed around the ai-axis (i.e., equal numbers oscillating in the a2- and a3-directions. 

Let P2, P3, and Pi be the power radiated by vibrational modes with electric dipole moment oscillating in the 
a3-, and ai-directions. For electric dipole radiation, we can imagine that each oscillator radiates 50% of its power in 
each of the two cardinal directions orthogonal to the dipole. 

For a nonrotating (or slowly rotating) grain, the angular momentum will undergo a random walk with 

where i, j, k run from 1 to 3 and the vector basis a^a^afe follow right-handed rule. 
Thus, defining J\ = Jf = J 2 , and noting that P2 and P3 are in-plane vibrations, and P\ is out-of-plane: 



d 72 , 1 



>Jt J±) ~ 2 



-N ip + N oop 



2h 2 , (j t Jf) = \n w 2K 2 (C23) 



D. COLLISION OF IONS WITH GRAINS 

The damping and excitation coefficients for the collision of ions with grain were calculated in DL98b, and refined 
in Ali-Hai'moud et al. (2009). The later study showed that the calculations by DL98b where the effects of dipole 
were disregarded remain valid. Thus, below we also disregard the effect of dipole in calculating the collision rate for 
ion-grain collision. The collision rate of ions with density rii with a grain of size a is 

f°° C°° { Tfl'V 2 \ 

Ri = ni J vTrbl^Anv 2 f(v)dv = m J virb^ATTv 2 Aexp \ ^-f J dv > ( D1 ) 

where A = (2feBT gas /mi7r) 3 ^ 2 , and 6 max is the critical impact factor so that all incoming ions with impact factor 
b < frmax collide with the grain. 

Dl. Neutral grain 

For neutral grain, & max is given by equation (|C17I) . thus, plugging (|C17I) into (|D1|) and performing the integral, we 
obtain 

1/2 

Ri(Z g = 0) = n^o? 



' 8k B T g3 

\ TO;7r 



1 , ^ ( 1Z 2 e 
2 V ak B T„, 



2„2 \ 1/2 



(D2) 



The rms angular momentum per collision is obtained by dividing the rms angular momentum to the collision rate. 
For the ion-neutral collision, we have 



/XT 2\ , rp 2 2 + 3 v / tt/2t + 2/t ak B T gas / a \f T gas \ 

(5 J ) = m k B T gas a - - , wfihr^— ~ = ™ (^ZT^J ) ■ 

D2. Charge grain 

For a charge grain, using the critical impact factor from (|C18[) for (|D1[) we obtain 

Rl (Z g *0) = n^a 2 (ZJ^A^g (Mf) , (D4 ) 
V m-xK J \ak B l gas J 

where g(x) = 1 — x for x < and g(x) = e~ x for x > 

The collision of one ion with the negatively-charged grain has rms angular momentum per collision for this case 
given by 

?- 9 2-2(b + (h 2 Z s e 2 , , 

(6 J 2 } = mi k B T gas a 2 V V , with <j> = % . D5 

1-0 ak B l„ as 
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Taking into account the grain charge distribution, the final form for collision rate becomes 

R ( 2Z 2 e 2 ^ 1/2 



R l = f(Z g = 0)n i na 2 



1 



2 V afcBTg as 



^ \ afc B Jg as ' 



Z s jtO 



where f(Z g ) is the grain charge distribution function. 

E. TRANSFORMATION OF COORDINATE SYSTEMS 

Damping coefficient A; = (AJj/Ai) and diffusion coefficients B.^ = (AJjAJj/Ai) are usually derived in the body 
coordinate system, while we are interested in the evolution of grain angular momentum in the inertial coordinate 
system. Let us define an inertial coordinate system eie2e3 in which the direction J is described by the angle £ between 
J with ei , and the azimuthal angle \ (see Figure 2 in Robergc & Lazarian 1999 in which angles </>, /3 are replaced by 
X, £) . To obtain these coefficients in the lab coordinate system, we first transform the body system to the external 
system xyz where z||J and x, y _L J (see Figure [TJ. Then, we perform the transformation from xyz system to the 
inertial system eie2e3. In the body system, the damping coefficients are given by 

* . <M) _ -Jl. _ JL ( ' ) , (ei) 

^£ ^gas,i 7"ed,z \ *Jjjn/B-^ gas / 

where i gas ,i = ^tot,i/TH,|| an d i = x,y and z with z\\&i. The diffusion coefficients, B b j = (AJ b AJ b /At) with = 
for i ^ j are related to the excitation coefficients as follows: 

B b = Bu = — — G tot ii, and B h = B b = Bj_ = — G to t,±- (E2) 

7h,|| yy T"H._L 

Transforming the vector A\ and matrix from the body system to the inertial system, and averaging over the fast 
precession of ai about J, we obtain damping coefficient (see Lazarian & Roberge 1997) 

Ai = <^7> = -— ^- { vYr ) , ** * = x,y,z , (E3) 

At Tgas.eff 7ed,eff \Oi||KB-i gas / 

where the effective gas damping reads 

T-gas.efiE = THJ-Ftotj '■ 2n ' ( E4 ) 

cos^t? + 7nsm 



and the effective time for electric dipole damping is given by 



TecLeff = T ed || — r -j-. (E5) 

cos 4 6> + 7edsm 9 

where 9 is the angle between ai and J, and 

7H = ^^, 7e d = ^ = / l 3 - (E6) 

If the coupling of parallel and per pen dicular rotation is accounted for, then the electric dipole damping is increased 
by a small amount, and equation (|E5j) is rewritten as (communication with C Hirata) 

Ted.eff - Tcd >H CQS 4 d + 7cdS i n ±0 + (ft3 + 3/l ) sm 2 6> COS 2 0/2" ^ 

We can see that for disk-like grains, th > 1, and r gaSj eff < th.|| ^tot,|| ? which results in faster gas damping time. 
The diffusion coefficients in the inertial system ej, B zz and B xx , B yy are given by 

B zz = B\\ Qsin 2 6>sin 2 £ + cos 2 (9cos 2 ^ + B ± Q [1 + cos 2 6>]sin 2 £ + sin 2 6»cos 2 ^ , (E8) 
B xx = B\\ Qsin 2 6»[cos 2 x + sin 2 xcos 2 £] + cos 2 6>sin 2 xsin 2 £^ + B ± Q[l + cos 2 6>] [cos 2 x + sin 2 xcos 2 £] + sin 2 6>sin 2 xsin 2 ^ , (E9) 

B yy = B\\ Qsin 2 0[sin 2 x + cos 2 xcos 2 £] + cos 2 6»sin 2 xsin 2 ^ + B± Q[l + cos 2 6»] [sin 2 x + cos 2 xcos 2 £] + sin 2 6>sin 2 xsin 2 ^ ,(E10) 

where £ is the angle between J and ei, and x is the azimuthal angle of J in the inertial system e^. 

In the presence of fast internal fluctuations, we need to average the damping and diffusion coefficients over 9. There- 
fore, the terms containing 9 in above equations are replaced by the averaged values, i.e., (cos 2 9) = J Q cos 2 8 Jlte {9) d9, 
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(sin 6>) = sin 2 6 f LT E(0)de. Note that we do not average the damping and diffusion coefficients over the precession 
angle x nere because the effect of magnetic field which results in the fast precession of J about the magnetic field is 
disregarded. 

Then, the Langevin equations in the inertial system become 

dJi = Aidt + \J Badqi, where (dq 2 ) = dt, and i = x, y, z. (Ell) 
F. DIPOLE EMISSION OF A TORQUE-FREELY ROTATING GRAIN 
Let us consider the simple case where the dipole is given by 

H = /Xyai (Fl) 
in the body system. In an external coordinate system xyz (see Figure [7]) , ai , a 2 and a3 are described as 

ai = sin sin 6*x — cos sin 9y + cos 9z, (F2) 
&2 = (cos <p cos ip — sin sin tp cos 6*)x + (sin cos ip + cos sin ip cos 9)y + sin ip sin 9z, (F3) 
a3 = — (cos sin ip + sin cos ip cos 6*)x + ( — sin sin ip + cos cos ip cos 9)y + cos ip sin 9z, (F4) 

where 0, ip and 9 are Euler angles (see Figure [7]). 

Let us now consider the torque-free motion of an axisymmetric grain with the angular momentum J and the ratio 
of inertia moments h. For this case, the symmetry axis ai precesses about J with constant angle 9 and the rate 0, 
and the grain rotates about the symmetry axis with the rate tp. They are given by (Landau & Lifshitz 1976) 

0=-=-, tp = —cos6(l-h). (F5) 
1± i\\ 

Precession and rotation of the grain with respect to J results in an instantaneous acceleration for dipole moment: 

fi = /iyai +/xxa 2 , (F6) 

where ai and a.2 are given by 

ai = — 2 sin 6>(sin 0x — cos 0y), (F7) 
— 2 (cos ^ cos "0 — sin (f> sin tp cos 9) — 2(pnp(— sin </> sin ip + cos <p cos tp cos 9) — "ip 1 (cos <p cos ip — sin cp sin tp cos 9 ) 

-0 2 (sin<?i>cos'0 + cos </> sin ip cos 0) — 20?/>(cos</!>sin'0 + sin cp cos ^ cos 9) — tp 2 (sin <p cos ip + cos (p sin ip cos 9) 
-ip 2 sinV>sin6'z, (F8) 
2 (cos (p sin -0 + sin cp cos tp cos 9) + 2<pip (sin <p cos tp + cos sin tp cos 9) + ip 2 (cos <p sin -0 + sin cp cos tp cos 9) 

<p 2 (sin cp sin ^ — cos <p cos cos 9) — 200(cos cos ip — sin sin ip cos 9) + ip 2 (sin sin ^ — cos cos ^ cos 9) y 
ip 2 cos ip sin 9z, (F9) 



a 2 = 



x 



y 



a 3 



Solving equation (|F5|) for Euler angles, and plugging into equation (|F6[) with the usage of equations (|F7|) and (|F9[) 
we obtain the dipole acceleration jl as functions of time. Performing Fourier transform for the components of fl gives 
us the spectrum and frequency of electric dipole emission. 

The dipole emission power of this torque-freely rotating grain can be obtained by averaging the fl 2 over Euler angles 
and ip in the range from to 2n: 

, s 2 . 2x f 2 * f 2 * d<P dip 2 , . 



REFERENCES 



Ali-Hai'moud, Y., Hirata, C. M., Dickinson, C, & Readhead. A. 

2009. MNRAS, 395, 1055 
Allers, K. N., Jaffe, D. T., Lacy, J. H., Draine, B. T., Richter, M. 

J. 2009, ApJ, 630, 368 
Anderson, N., & Watson, W.D. 1993, A&A, 270, 477 
Bennett, C. L. et al. 2003, ApJS, 148, 97 

Bouchet, F. R., Prunet, S.,& Sethi, Shiv K. 1999, MNRAS, 302, 
663 

Boughn, S. P., & Pober, J. C. 2007, ApJ, 661, 938 
de Oliveira-Costa et al. 1999, ApJ, 527, 9 
de Oliveira-Costa et al. 2002, ApJS, 567, 363 

Dobler, G., Draine, B., & Finkbeiner, D. P. 2009, ApJ, 699, 1374 



Draine, B.T., & Lazarian, A. 1998, ApJ, 494, L19 (DL98a) 
Draine, B. T., & Lazarian, A. 1998, ApJ, 508, 157 (DL98b) 
Draine, B. T., & Li, A. 2001, ApJ, 551, 807 
— . 2007, ApJ, 657, 810 (DL07) 
Draine, B.T., & Sutin, B. 1987, ApJ, 320, 803 
Efstathiou, G. 2003, MNRAS, 346,26 
Erickson, W,C, 1957, ApJ, 126, 480 
Ferrara, A., & Dettmar, R.-J. 1994, ApJ, 427, 155 
Finkbeiner, D. P., Langston, G. I., & Minter, A. H. 2004, ApJ, 
617, 350 

Gardiner, C. W. 1983, Handbook of Stochastic Method (Berlin: 
Springer) 



27 



Gold, B., Bennett, C. L., Hill, R. S., et al. 2009, ApJS, 180, 265 
Gold, B., Odegard, N., Weiland, J. L., et al. 2010, ApJS, 
submitted 

Hoang, T., & Lazarian, A. 2008, MNRAS, 388, 117 
Hoang, T., & Lazarian, A. 2009, ApJ, 695, 1457 
Jones, R.V., & Spitzer, L. 1967, ApJ, 147, 943 
Kogut, A. et al. 1996a, ApJ, 460,1 
Kogut, A. et al. 1996a, ApJL, 464, 5 

Landau L.D, & Lifshitz E. M. 1976, Mechanics (Oxford: 

Perganon) 
Lazarian, A. 1994, MNRAS, 268, 713 
Lazarian, A. 1995, MNRAS, 274, 679 

Lazarian, A. 2007, J. Quant. Spectrosc. Rad. Trans., 106, 225 
Lazarian, A., & Draine, B.T. 1999, ApJ, 520, L67 
Lazarian, A., & Hoang, T. 2009, arXiv 0901.0146 



Lazarian, A., & Roberge, W. 1997, ApJ, 484, 230 
Li, A., & Draine, B. T. 2001, ApJ, 554, 778 

Mathis, J.S., Mezger, P.G., & Panagia, N. 1983, A&A, 128, 212 
Purcell, E.M. 1979, ApJ, 231, 404. 

Ragot, B.R. 2002, ApJ, 568, 232 Roberge & Lazarian 1999 

Roberge, W., DeGraff, T.A., Flaherty, J.E. 1993, ApJ, 418, 287 

Silsbee, K., Ali-Haimoud, Y., & Hirata, C, arXiv 1003.4732vl 

Sironi, L., & Draine, B.T. 2009, ApJ, 698, 1292 

Tegmark et al. 2000, ApJ, 530, 133 

Ysard, N.,Verstraete, L. 2009, arXiv: 09063102 

Weingartncr, J. 2009, ApJ, 690, 875 

Wcingartner, J.C., & Draine, B.T. 2001b, ApJS, 134, 263 



